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SYMBOLS 


a 


half length of notch 


c 


superscript indicating compressive failure stress 


tensile modulus of laminate 


^ll'^22 


tensile moduli of lamina in fiber and transverse 
directions , respectively 


xy 


axial shear modulus of laminate in xy coordinate 
system 


’x 'y ' 


axial shear modulus of laminate in x'y' off-axis 
coordinate system 


’12 


axial shear modulus of lamina 


H 


laminate thickness 


h 


lamina thickness 


k 

m 

N 


n 


- empirical constant used in Equations (3) and (4) 

- size of overstressed intact region 

- number of cycles 

- fatigue lifetime 

- superscript indicating notch stress 


R 


a . /a , stress ratio 
min' max ' 


S 


cyclic maximum stress/static failure strength 



- cyclic maximum shear stress/static failure shear 
strength 

- effective combined stress state 

- superscript indicating tensile failure stress 




non-dimensional and dimensional axial displacement 
inside the core region, respectively 


Ui,Ui 


non-dimensional and dimensional axial displacement 
in the over stressed region, respectively 


'*2 '^2 ” non-dimensional and dimensional axial displacement 

in the average uniformly stressed region, respec- 
tively 


V. 



X 


direction of load in laminate 


y 


transverse direction to the load 


X 


,Y' 


direction of and transverse direction to an 
off-axis crack in a laminate 


a 


Y 


ult 


- non-dimensional axial "crack" length (including 
plastic zone) , inelastic length 

- axial shear failure strain 



a 

max 

°min 

'^r 

°SCFM 



- non-dimensional axial crack length 

- lamina orientation 

- applied gross laminate stress in x direction 

- applied laminate stress at which axial crack 
propagates to infinity 

- initial static 

- axial tensile strength in x'y' coordinate system 

- applied laminate stress at which axial crack 
is initiated at notch tip 

- applied laminate stress when material adjacent to 
the notch reaches 

SCFM 

- maximum tensile cyclic stress 

- minimum tensile cyclic stress 

- residual strength 

- maximum composite overstress in material adjacent 
to notch 

- applied laminate stress at which transverse crack 
propagates from notch tip 

- initial static 

- laminate axial shear stress 

- unnotched laminate tensile strength 


unnotched laminate compressive strength 
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^ >1 


laminate stress components in x*y* coordinate 
system 


0 

X 


I 9 


o 
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I 9 


x'y' 



0 


n 




T ' 

max 



I 


T 


n 


T 


n* 


n 

e 

d 


T® 

amax 


applied stress at which axial plasticity is 
initiated at the notch tip 

compressive and tensile strengths of the unnotched 
laminate in y* direction, respectively 

1/2 (a°, - Oy.) and 1/2(0°, + o^,) , respectively 

notch root tensile stress in xy coordinate system 

due to CT_.^ 
max 

stress components in lamina 

laminate strength for off-axis crack failure mode 

maximum cyclic shear stress level 

shear yield. stress and strength of the unnotched 
laminate in the xy coordinate system 

layer axial static shear strength 

shear yield stress and strength of the unnotched 
laminate in the x'y' coordinate system 

effective laminate shear yield stress and strength 
due to the combined stress state in the x'y' co- 
ordinate system 

notch root shear stress in xy coordinate system 
due to 

max 

notch root shear stress in x'y' coordinate 
system due to only 

effective laminate notch shear stress 


reduced shear strength at notch root used in 
Eq. (1) 

increased notch shear stress used in Eq. (5) 

laminate axial shear strength which correlates 
with axial crack propagation failure at a in 
a notched laminate 
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FATIGUE OF NOTCHED FIBER COMPOSITE LAMINATES 
PART I: ANALYTICAL MODEL 

By P. V, McLaughlin, Jr., S. V. Kulkarni, 

S. N. Huang and B. Walter Rosen 
Materials Sciences Corporation 

SUMMARY 

This report describes a semi-empirical, deterministic analy-. 
sis for prediction and correlation of fatigue crack growth, resid- 
ual strength, and fatigue lifetime for fiber composite laminates 
containing notches (holes) . The failure model used for the 
analysis is based upon composite heterogeneous behavior and exper- 
imentally observed failure modes under both static and fatigue 
loading. The analysis is consistent with the wearout philosophy. 

During fatigue loading of a notched laminate, high axial 
shear stresses exist near the notch root. These high shear 
stresses will cause the material in the vicinity of the notch to 
be degraded much more rapidly than throughout the rest of the 
Iciminate. As the properties of laminae change in the notch region 
due to cyclic loading, the strength of the notched laminate will 
alter. Residual strength of a notched composite is therefore 
treated as a static failure utilizing fatigue-degraded lamina 
properties in the notch region. A fatigue failure occurs when 
the residual strength of the notched composite falls to the level 
of applied maximum cyclic stress. 

Laminate axial shear behavior governs notched composite 
fatigue behavior. The fatigue analysis method therefore consists 
of the following: 

(1) Estimating the notch root stress state for a given 
laminate tensile stress. 

(2) Performing a laminate analysis to obtain lamina stresses. 

(3) Computing fatigue-induced material property changes in 
each layer . 



(4) Predicting new laminate properties in axial shear. 

(5) Calculating the changed residual strength properties of 
the notched composite. 

The process is repeated as necessary utilizing convenient incre- 
ments of number of cycles. 

The modes of failure treated by the present analysis are: 

(1) Axial cracking in the load direction. 

(2) Transverse cracking across the specimen. 

(3) Cracking at an angle to the load axis along a fiber 
direction. 

Axial cracking and transverse cracking failure modes are treated 
together in the analysis. Cracking off-axis is handled by making 
a modification to the axial cracking analysis. 

The analysis predicts notched laminate failure from unidirec- 
tional material fatigue properties using constant strain laminate 
analysis techniques. For multidirectional laminates, it is nec- 
essary to know lamina fatigue behavior under axial normal stress, 
transverse normal stress and axial shear stress.' 

Examples of the analysis method are given for [0] and 
[ 02 / 145 ] boron/epoxy laminates with varying notch sizes. Results 
show that static failure and fatigue failure modes need not be 
the same. This fact suggests that the often-made assumption of 
one-to-one correlation between static and fatigue failure distrib- 
utions should be critically reexamined. 

The present model has not been experimentally verified for 
fatigue loading, although good correlation has been obtained for 
static failure. Also, the present model does not include the ef- 
fects of stacking sequence or interlaminar debonding. Further 
work should concentrate in these areas . 



INTRODUCTION 


The utilization of any material for primary aerospace 
structural applications requires adequate safety margins under 
all anticipated environments and an adequate structural lifetime. 
This requires a sound understanding of the fracture and fatigue 
behavior of the structural material. Over the years, a great 
body of empirical data and analytical techniques have been gen- 
erated for metallic alloys. In contrast, relatively little is 
known about the fatigue behavior and residual strengths of compos- 
ite materials, particularly when the material has a notch or 
other source of stress concentration. Furthermore, for any 
composite material (e.g., boron/epoxy) , there are many aspects of 
the internal material geometry which can be varied, such as 
amounts of fibers oriented in different directions. Thus, the 
development of a broad experimental data base to completely 
characterize one material system is a costly undertaking. For 
this reason, an analytical model of the fatigue behavior of com- 
posite laminates is desirable for the design engineer. 

Until recently, fatigue analysis took the form of empirical 
generation of standard S - N curves and perhaps residual stress 
information for specific materials and particular composite lay- 
ups. An important addition to fatigue thinking was made with the 
conception of a "wearout" philosophy of composite material degrada- 
tion under cyclic loading (Reference 1) . This philosophy treats 
fatigue damage as the growth of preexisting flaws or discontin- 
uities in a material. As these flaws are made larger due to the 
action of cyclic stresses, there is a generally greater tendency 
to failure and a lower residual strength. When residual strength 
drops to the level of maximiam cyclic stress, fatigue failure 
ensues. If a model can be developed for the growth of these flaws 
and their effect upon static strength, then residual strength and 
lifetime of materials can be predicted. 

This report presents a semi-empirical analysis for prediction 
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and correlation of fatigue crack growth, residual strength, and 
fatigue lifetime of fiber composite laminates containing notches. 
The analysis utilizes a material model which treats physically 
realistic failure modes. The approach is consistent with the 
wearout fatigue philosophy, in that the layered material prop- 
erties degrade with cyclic loading and cause reduced residual 
strength and eventual fatigue failure of a notched composite mate- 
rial. A major difference between the current analysis and that 
of Reference 1 is the mechanics of crack propagation employed. 
Reference 1 treats crack propagation in composite materials using 
concepts of classical fracture mechanics which have been developed 
for homogeneous materials. The present work, on the other hand, 
incorporates the main physical characteristics and failure mech- 
anisms of fiber composite materials which have been observed 
during static and fatigue tests (see, for example. Reference 2) . 
This topic is given further discussion in the next section. 

The philosophy behind the present fatigue analysis is that 
cyclic loading of a notched composite causes localized property 
degradation in the vicinity of the notch. This combination of 
reduced properties and stress concentrations due to static loading 
causes localized damage in the vicinity of the notch. This 
damage results in a redistribution of local stresses and hence 
a change in strength. When these fatigue-induced changes cause 
a decrease in residual strength with number of cycles N, the 
fatigue lifetime is the number of cycles at which the resid- 
ual strength has reduced to the applied stress. Implicit in 
this philosophy is the notion that a static failure model will 
always be able to determine strength and stress-strain behavior 
if the instantaneous properties of the material are known at any 
given N. Therefore, the two main ingredients in this fatigue 
analysis model are a valid static failure model plus a method 
for determining the degradation of material properties and growth 
of flaws with number of cycles. The approaches taken to these 
two topics under the present research program are described 
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following a discussion of analytical and experimental considera- 
tipns essential to the fatigue analysis model. The body of this 
report contains material which is judged essential to under- 
standing the analysis. Mathematical details are presented in 
appendices . 
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ANALYTICAL AND EXPERIMENTAL OBSERVATIONS 


The behavior of notched fiber composite laminates under 
static and fatigue loading exhibits many features which are not 
observed for homogeneous materials such as metals (Ref. 2, 5, 6). 
For example: 

(i) Frequently, the direction of crack growth in a fiber 

composite laminate is not perpendicular to the applied 
load, even for perpendicular slit notches. 

(ii) Static failure modes may be completely different from 
the mode of failure observed during or after a program 
of cyclic loading. 

(iii) The residual tensile static strength of certain notched 
composite laminates subjected to fatigue loading can be 
greater than the initial static strength of those lam- 
inates . 

Due to these observed inconsistencies with homogeneous mate- 
rial fatigue and fracture mechanics theories, it is not expected 
that classical theories will have general application to composite 
laminates on the macroscale. Limited success has been obtained 
in applying classical fracture mechanics to composites when the 
mode of failure is a collinear crack growth from a preexisting 
notch (Reference 3) . Unfortunately, composites exhibit several 
types of non-collinear crack growth which cannot be treated by 
classical means. Several attempts have been made to utilize 
classical fracture mechanics by treating fiber and matrix as 
distinct homogeneous materials and analyzing details of stress 
and strain around cracks or flaws in the separate phases (e.g.. 
Reference 4) . While these studies provide insight into the 
mechanisms of crack growth in heterogeneous materials, they are 
not well developed and appear to be too complicated for inclusion 
in an overall analysis of composite laminate fatigue behavior. 

Tests of many fiber composite material systems (e.g., see 
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References 2 and 5) have shown that there are primarily three 
modes of failure in a notched fiber composite laminate under ten- 
sion. All three failure modes have been experimentally observed 
in both static and fatigue tests. The failure modes are illus- 
trated in Figure 1 which contains photographs from Reference 2 
of failed boron/epoxy laminates that have undergone fatigue cycling. 
An axial crack propagation mode (Figure la) occurs when the central 
portion of the laminate near the notch pulls out of the remaining 
material due to high shear stresses in the notch vicinity. The 
axial crack propagates to the grips, ultimately causing specimen 
failure. A transverse or collinear crack propagation mode is 
shown in Figure lb. Here, material immediately adjacent to the 
notch has been over stressed to the' point where an unstable crack 
continues to grow perpendicular to the load direction. If there 
are fiber families oriented at an angle to the direction of load- 
ing (Figure Ic) , a crack may initiate and extend parallel to 
these off-axis fibers, resulting in an off-axis crack propagation 
failure mode. As shown in Figure Id, failure can occur by sev- 
eral of these modes in the same laminate. 

In all failures, it has been observed that cracks propagate 
either along a fiber direction or transversely across the specimen 
perpendicular to the direction of tensile loading. The mode of 
failure can be different for fatigue or static loading of the same 
laminate. In addition, different notch sizes can trigger dif*- 
ferent failure modes in the same laminate. Other factors which 
can affect failure mode are ply orientation and constituent mate- 
rial properties. 

Any fatigue failure model of a fiber composite laminate must , 
exhibit the above described failure modes , preferably through a 
built-in recognition of the key aspects of fiber composite mate- 
rial heterogeneity. The model described in this report has 
attempted to do this in a way which is sufficiently general to 
apply to all fiber composite laminates. 
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STATIC FAILURE MODEL 


The primary modes of both static and fatigue failure of 
notched composites are; 

(a) Propagation of a crack parallel to one of the fiber 
directions (either axial or off-axis) . 

(b) Propagation of a crack transversely across the specimen 
normal to the direction of applied tensile stress. 

The basic static model for axial and transverse failure of notched 
composites was developed in Reference 6. Modifications to the 
basic model and addition of off-axis cracking was performed 
during the present program. Because of its importance to the 
fatigue analysis, a discussion of the static failure model and 
examples for all modes of failure will be given here. 

The basic model for static failure of notched composites 
treats axial and transverse crack propagation. Axial failure is 
caused by the tendency of the material separated by the notch to 
"pull out" . This causes high shear stresses in the direction of 
load and may result in propagation of an axial crack. Transverse 
crack propagation is caused by the stress concentration effect 
of the notch. Material adjacent to the notch may be overstressed 
to failure, causing a crack to run perpendicular to the load. 

When a crack propagates at an angle between 0® and 90° to 
the load direction, it does so under the action of combined 
shear and normal stresses along the crack. The off-axis phenom- 
enon has been judged to be physically similar to axial crack 
propagation, but with the added complexity of a combined stress 
state. The off-axis crack propagation model has, therefore, been 
developed from the axial crack propagation model by the inclusion 
of combined stress effects. 

The analysis uses the same model to treat axial and transverse 
crack propagation, and a modified model to treat off-axis crack 
propagation. In the following subsections, the axial/transverse 
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failure model and the off-axis failure modification are summarized. 
Mathematical details can be found in Appendix D. 

Axial and Transverse Crack Propagation Model 

The model for static failure of a notched composite laminate 
is illustrated in Figure 2 , and coordinate system conventions are 
shown in Figure 3. The laminate is asstamed (Figure 2) to be 
under a tensile load a in the x direction. A notch of width 2a 

X 

is centered in the specimen at x = 0. The coordinate y lies in 
the laminate plane transverse to the direction of applied tensile 
load. Reference 2 has shown that notch shape is not as important 
to failure as the size (width) of the notch. Hence, a slit, cir- 
cular hole, or square hole of width 2a will all be modeled in a 
y direction discontinuity in the laminate of 2a. Unless specif- 
ically noted to the contrary, the word "notch" will apply to any 
through-the-thickness hole of any shape. 

In the central core region (Figure 2) , the core wants to 
"pull out" from the notch area due to the applied tensile loading. 
This core, however, is restrained by shear stresses between the 
core and the adjacent intact material. These shear stresses 
generally result in a region of high shear strain parallel to the 
loading direction along the 0° fibers. Immediately adjacent to 
the notch core, there is a region of width m which is overstressed 
due to presence of the notch. It is assumed that this region 
bears the full stress concentration effects. Everywhere outside 
the core and overstressed regions, the laminate is uniformly 
strained. This uniformly strained region is called "average 
material", because there are no overall average stress or strain 
gradients and the laminate is assumed to behave as a homogeneous 
anisotropic material. 

Shear strain due to core pullout is assumed to be effective 
over a region three times the size of the overstressed fiber 
regions. This shear strain region extends a distance m into the 
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region of uniformly strained average material which lies adjacent 
to the overstressed region. The total thickness of the laminate 
is H, and the thickness of each layer in the laminate is h. 

Due to the high axial shear strain region near the notch, 
description of the behavior of the laminate in axial shear, 
is very important. The shear stress-strain curve can vary any- 
where from nearly linearly elastic, as would occur for a 
[0 . /±45 ./90, ] laminate to one which is highly nonlinear, such as 

3. 3 X 

would occur for a [0] laminate. The axial shear stress-strain 
curve is therefore approximated as linear elastic-perf ectly 
plastic (Figure 4) . The shear yield stress is assumed equal to 
the shear strength, t°. The axial shear modulus, is chosen 

to obtain the best fit to the actual laminate shear stress- 
strain curve. 

As the gross laminate tensile stress a is increased, the 
shear stresses at the notch tip will increase. When the yield 
stress in shear t° is reached at the notch tip, continued loading 
will cause an inelastic region to form in the axial direction. As 
the axial shear strain at the notch tip continues to increase with 
increased stress, the axial shear failure strain will be 

exceeded and a crack will begin to grow in the axial or x direc- 
tion. When this crack grows very large (mathematically, to infin- 
ity) the notched laminate is considered to have failed in an axial 
crack propagation mode. If, however, material in the overstressed 
region adjacent to the notch is stressed past the tensile strength 
of the unnotched laminate, a crack will propagate in the transverse 
or y direction. In this latter case, the laminate is considered 
to have failed in a transverse or collinear crack propagation 
mode. 

Appendix A fully describes the static failure model and re- 
sulting stress and deformation analysis for axial and transverse 
crack propagation modes. It is, in essence, an inelastic shear 
lag analysis with the axial shear behavior of the laminate playing 
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9-n extremely important role. 

In the analysis, ? is a nondimensional axial crack half- 
length extending from the notch tip (x = 0) to the bottom of the 
plastic zone. The non-dimensional distance from the notch to 
the beginning of the elastic zone is a, and includes axial crack 
and plastic zones (see Appendix A for details of the non-dimen- 
sionalization) . For this reason, a is called the "inelastic 
length" and is a good measure of the extent of damage in the 
axial direction from the notch tip. 

The analysis has the capability of computing both the in- 
elastic length and maximim overstress in material adjacent to the 
notch for a given applied laminate tensile stress so that tendency 
to axial or transverse crack propagation modes can be monitored. 
Typical results of the analysis and a demonstration of the proced- 
ure to be followed in computing axial and transverse failure 
stresses .are presented below. 

Analysis procedure and typical results fCr axial and transverse 

cracking .- Typical results of the static failure analysis for a 

notched laminate are shown in Figure 5. Figure 5a illustrates the 

growth of the damaged region with increased applied tensile stress a. 

Until the yield stress of the laminate in shear is exceeded at the 

notch tip, a will be zero. When the yield shear stress t° is 

exceeded at the notch tip, the inelastic length a will begin to 

grow. The .applied stress at which plasticity is initiated at the 

notch tip is called a . At some higher stress level a the shear 
^ y c 

strain in the laminate at the notch tip will exceed 
laminate shear failure strain. An axial crack will then begin to 
grow from the notch tip. As loading progresses, the damage zone 
will increase without bound as the applied stress assymptotically 
approaches a level a^. This stress represents the failure stress 
of the notched composite in the axial crack propagation failure 
mode . 
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As stressing in the composite increases and as the inelastic 
length a grows, the region adjacent to the axial crack is over- 
stressed by varying degrees (Figure 5b) . As the applied stress 
is increased, it is possible that the maximum overstress in the 
adjacent material, o , will increase to the point where it 
exceeds the unnotched laminate strength, In this case, 

transverse crack propagation ensues. The applied laminate stress 
at which transverse cracking occurs is called it is also 

possible that, as the axial crack grows, a blunting of the stress 
concentration effects of the notch will occur. This could cause 
a decrease in stress in the adjacent material before it reaches 
failure levels, and remove any possibility of transverse crack 
propagation occurring. 

The analysis is capable of computing both and assuming 
that neither failure mode has yet occurred. If the computed 
value of CJ^ is lower than a^, axial cracking is the predicted 
mode of failure. In this case, is the notched composite 
failure stress. If, on the other hand, the computed cr^ is lower 
than a , transverse cracking is the predicted mode of failure. 

A 

The failure stress of the notched composite in this case is a^. 

Discussion .- Two key ingredients to the failure predictions 
are the static stress-strain behavior in axial shear of the lam- 
inate and the unnotched axial tensile behavior of the laminate. 

These quantities can either be determined from laminate experimental 
data or from state-of-the-art laminate analysis techniques provided 
that layer data is available. In addition, the size of the shear 
transfer zone next to the notch, 3m, and the size of the stress 
concentration region, m, are of prime importance to the analysis. 

The basis for their selection is discussed in the following par- 
agraphs . 

The size of the shear transfer zone adjacent to the stress 
concentration region was chosen to be three times the size of the 
stress concentration region based on physical reasoning. There 



the axial displacements inside the core region (U ^) , in the 
overstressed region (U^^) , and in the flanking uniformly 
strained material (U 2 = a^re generally different from one 

another (see Figure 2) . Shear stresses will therefore be develop- 
ed between the core region and stress concentration region, and 
between the stress concentration region and the uniformly stressed 
regions. If the differences in displacements and 

are of the same order of magnitude, it is likely that shear . 
strains will extend an equal distance into the core and into the 
average material. The choice of stress concentration region size 
m for these shear strain penetration distances was based on the 
stress concentration region size itself being a measure of stress 
perturbations near the notch. 

The size of the overstressed region, m, is a parameter which 
must be experimentally determined for a given laminate. It is 
calculated from the static tensile strength of the laminate having 
a notch of a given size. The parameter m is chosen such that the 
analytical prediction for notched laminate strength and the mode 
of failure thereof match experimental data. Once the value m has 
been determined for one notch size in a particular laminate, it 

1/2 

can be computed for any notch size from the equation m/m<,= (a/ao) 
where mo is the experimentally determined value of m from tensile 
strength data on a laminate with notch width 2ao . This square 
root relationship is an experimentally determined correlation from 
tests on boron/epoxy laminates, and is described in Reference 6. 

Off-Axis Crack Propagation Model 

The preceding subsection has treated failure of notched lam- 
inates when the mode of cracking is either axial along the 0® 
fiber direction, or transverse to the specimen normal to the 
direction of loading. For some Iciminates, such as a [±45] boron/ 
epoxy laminate, the mode of fatigue cracking is off-axis along a 
fiber direction (in this case, the ±45 degree direction) . For 
this mode of failure, the preceding model and analysis are not 
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directly applicable without modification. This subsection 
describes modifications to the original model which were made to 
allow analysis of off-axis cracking . 

Notch inclination and combined stress effects .- Figure 6 
illustrates the stress state around a notch which is tending to 
develop a crack along the ±0 fibers. Along the load axis the gross 
laminate stress state is axial tension. However, the crack is 
not growing in the direction of load. Suppose material axes 
x' and y' are selected in the 0 direction (Figure 6). The crack 
will be growing in the direction of the x' axis tensile stress 
a I. Due to the rotation of axes, the tensile stress a , will 
not be the same as the x direction stress a. In addition, there 
will generally be substantial axial shear and transverse tensile 
stresses in the x'y' coordinate system. Also, the notch is now 
inclined at an angle 0 to its original position. The following 
paragraphs describe how notch inclination and combined stress 
states are analyzed so that off-axis cracking can be handled 
directly by the original model referred to the x'y' coordinate 
system. 

It has been shown (Reference 2) that the shape of the notch 
(slit, rectangular, circular, etc.) in a laminate does not 
greatly affect composite failure. The width of a notch perpendic- 
ular to the direction of loading is the only notch dimension 
which is of importance. Therefore, the slit notch configuration 
shown in the xy and x'y' coordinate systems (Figures 6a and b, 
respectively) , can be treated for analysis purposes as shown in 
Figure 6c. In this case, the original notch has been replaced by 
an equivalent notch of size 2a cos 0 in the y' direction. (Note 
that no change need be made for a circular hole, and that an 
increase or decrease in equivalent notch size may be necessary 
for a rectangular hole.) 

Performing an exact analysis of axial cracking in the presence 
of a combined stress state is extremely difficult. Therefore, 
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combined stress effects are taken into account in an approximate 
fashion. Laminate gross transverse tension Oy, and axial shear 

are superimp, osed on the axial shear stress at the notch root 
which would be computed from the axial crack analysis due to only 
(see Figure 7) . The main effect of this superposition will be 
to increase the likelihood of x'y' axial shear failure. Through- 
out the loading process, the axial shear stress reduces the 

stress required from the shear-lag analysis to cause axial cracking. 
The transverse normal stress , is assumed to further reduce the 
shear strength. 

The assumed effect of axial shear and transverse stresses is, 

therefore, to change the laminate axial shear strength inputted 

into the standard tensile strength analysis. One method of doing 

this is presented in Appendix B. The unnotched laminate shear 

strength t° ' in the x'y' coordinate system is assumed to be directly 

reduced by applied shear a , , . A quadratic interaction is chosen 

X y 

for the reduction in x° ' due to result .is a reduced, or 

effective, unnotched laminate shear strength ' which depends 
upon the level of applied stress a , through cr , , and o Fig- 

ure 8 shows a typical variation in t_® ' with a , . 

6 X. 

Procedure for off-axis failure analysis .- If the effective 
shear strength were a constant, it could be used directly in the 
failure model to compute r the value of o^, which causes failure. 
The applied failure stress Cg required for off-axis failure could 
then be obtained from stress transformation equations with a^, = 
and the corresponding values of a^., and Oj^iyi at failure. However, 
Tg“ ' varies with applied load and a different procedure must be 
adopted. 

First the x' crack failure stress a^i is computed from the 
original model for a range of effective shear strength Xg® ' which 
goes from near zero to x° ' (the x' direction unnotched laminate 
shear strength). From this data, a curve of vs. Tg® ' can be 
constructed as shown in Figure 8. Next, the variation in Xg®' as 
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a function of a^, is computed from the Appendix B equations. This 
information is plotted on the same graph as the a,, vs. x curve 
in Figure 8 with and a^, axes coincident. The intersection of 

the two curves defines both ' and the failure stress . 

Once a^, has been determined, the corresponding stress state 
failure is known. The stress transformation 

X y * y 

equations are then used to find the corresponding tensile failure 
stress Oq in the original xy coordinate system. 

Static Failure Analysis Procedure 

The static failure analysis procedure for notched composites 
which has been developed utilizing the axial/transverse model and 
the off-axis model is as follows: 

1. The axial/transverse failure model is exercised to compute 
the laminate stresses causing axial crack failure (a^) 

and transverse crack failrue (a^) . 

2 . The off-axis model is exercised to compute the laminate 
stress Og which causes failure due to a crack running 
along 6-direction fibers. 

3. The failure stresses cr^, and Oq are compared. The 

lowest is the predicted failure stress and represents the 
dominant failure mode. 

Examples of the static failure analysis procedure are given 
in Appendix C for boron/epoxy unidirectional and multidirectional 
laminates. Figure 9 presents the procedure in more detail. 

Discussion 

It is noted that axial and off-axis crack propagation failure 
stresses may not be "ultimate" stress predictions in the usual 
sense for a uniaxial tensile specimen. This is especially true if 
there are elaborate load transfer devices at the grip. Once axial 
cracks have propagated to the grips during a tensile test, the 
specimen has become two smaller unnotched specimens. If the axial 
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crack is clean, loading may continue to be applied, and the spec- 
imen may fail at a load equal to the unnotched laminate strength 
times the net cross-sectional area. In a structure, however, large 
scale axial crack propagation would be considerably more serious 
because of the lack of crack arresting mechanisms such as the tabs 
on a tensile specimen and the grips of the testing machine. In 
addition, structural panels seldom, if ever, experience only 
uniaxial tensile loading. Under the presence of even small combined 
stresses, a large axial crack will precipitate catastrophic failure. 
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FATIGUE MODEL 


Fatigue failure in notched fiber composite laminates occurs 
as a result of growth of cracks along preferred directions in the 
laminate as load cycling proceeds. This crack growth is influenced 
by local material properties which change as a result of the re- 
peated loads. When the local material properties and the geometry 
of the damaged regions are defined after any number of cycles, the 
residual static strength of the material can be determined from the 
static model outlined in the preceding section. Thus, as a first 
requirement, it is necessary to determine the material properties 
of the laminate in the notch vicinity as functions of the niimber 
of load cycles. With this information, the growth of the damaged 
region can be defined. 

The main modes of damage growth as discussed previously are: 
crack propagation along the load axis , crack propagation trans- 
verse to the loading direction, and crack propagation off-axis 
along a fiber direction. All of these crack and failure mech- 
anisms are determined on a static basis by the analysis presented 
in the preceding chapter. The main inputs to the static failure 
analysis are the laminate elastic properties and the laminate 
failure stresses in axial tension and axial shear. As these 
properties change with cyclic loading, the static failure behavior 
of the notched composite will alter. If the properties are known 
at any given cycle, it is possible to perform a static analysis to 
determine the static residual strength of the notched composite. 

If the residual strength so computed is less than the maximum 
cyclic stress, the fatigue lifetime of the notched composite has 
been exceeded. 

The main ingredients to the fatigue analysis, then, are: 

(a) Determination of changes in material properties in the 
notch region with number of cycles, and 

(b) Utilization of these properties in a static failure 
analysis to predict residual strength. 
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The approaches taken to these aspects are discussed in the fol- 
lowing subsections . 

Notch Region Material Degradation 

In a notched composite, the main cause of crack growth is the 

tendency of the notched core region to pull out, generating stress 

concentrations in adjacent material and very high shear stresses 

in the axial direction. The remainder of the composite outside 

the notch region is generally not highly overstressed. Therefore, 

most fatigue degradation will occur due to shear and tensile 

stresses in the region immediately surrounding the notch, and the 

main result will be to alter the axial tensile and axial shear 

behavior. Specifically, it is anticipated that the axial Young's 

modulus, E ,:the unnotched tensile- strength,, g^, the’ shear 

modulus G , the axial shear failure strength, t“, and the laminate 
xy 

axial shear failure, strain will be primary laminate properties 

which vary under cyclic loading. Due to the high shear stresses 
near the notch, the shear properties will degrade at a higher rate 
than the tensile properties. An assumption to this effect is not 
necessary, however, as will be seen later. 

The most direct way to determine laminate tensile and shear 

behavior under cyclic loads is simply to perform fatigue tests on 

unnotched specimens. This method, though straightforward, would 

# 

require a battery of tests for each laminate layup geometry in 
both the axial xy and off-axis x'y' systems. The results would 
only be applicable to the specific laminate tested; each new lam- 
inate would require its own battery of tests. 

Clearly, if a method of generating Icuninate fatigue behavior 
from lamina properties can be derived, the number of basic property 
tests required can be reduced. Lamina properties so-generated 
would be applicable to laminates with any layup geometry and loaded 
in any coordinate direction. The method chosen for this analysis 
utilizes constant strain laminate analysis techniques, and is as 
follows : 
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1. The laminate stress state which exists in the notch vicin- 
ity (combined axial tension and axial shear x^) is 
determined. The method for so-doing is outlined later and 
presented in detail in Appendix D. 

2. A constant strain laminate analysis is performed on the 

laminate under the stresses a and x utilizing initial 
static lamina properties. From this laminate analysis, the 
axial normal stress, ^^. 1 ' axial shear stress <^2.2' 

the transverse normal stress, 0^2 which exist in each . 
layer of the laminate are computed. These stresses 
represent the maximum stress due to cyclic loading in 
each lamina. 

3. From data on fatigue behavior of a unidirectional lamina, 
calculations are made of changes in lamina elastic prop- 
erties and lamina residual failure stresses after some 
small increment in cycles. It is assumed that, over the 
range of cycles considered, the individual stresses in 
each lamina are not significantly changing. Therefore, 
the degradation which occurs in that cyclic loading 
interval does so under a constant stress state. This is 
an approximation, as it is expected that small changes 
will occur in each cycle. However, the error in such a 
"step function" approximation can be minimized by choosing 
very small increments of N for computation. 

4. The changed lamina elastic properties and lamina strengths 
are re-introduced into the laminate analysis to predict 
new laminate elastic properties and failure stresses. 

5. A second increment of cyclic loading is selected, and 
lamina stresses which are occurring in the notch vicinity 
are re-computed using the new lamina and laminate prop- 
erties. The resulting lamina stresses are used with 
unidirectional material fatigue behavior to predict 
further changes in lamina properties which may have 
occurred over the second increment of cyclic loading. 

The new lamina properties from the second increment of 
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cyclic loading are then used in the laminate analysis to 
generate laminate elastic properties and residual failure 
strengths in the notch vicinity. 

This procedure can be repeated enough times as is necessary 
to find the residual strength after a required number of cycles, 
or until the residual strength of the laminate reduces to the 
stress levels existing at the notch. When the latter occurs, the 
fatigue lifetime of the material near the notch has been reached. 

Figure 10 illustrates the kind of lamina fatigue information 
which would be necessary to generate laminate fatigue properties 
in the notch vicinity. Only axial tensile strength, axial shear 
strength, and axial shear modulus are illustrated, although 
transverse Young's modulus information may also be necessary. 

This topic is given further discussion in ' later subsections of 
this report. 

Figure 11 illustrates typical notch root material fatigue 
information which would be generated from lamina fatigue behavior 
using the above laminate analysis approach. Only axial shear 
information is shown, but axial tensile properties can also be 
computed. The following subsections outline the use of this 
information with the static failure model to perform fatigue analy 
sis of notched composite laminates. The axial/transverse fatigue 
modes are discussed first, and the off-axis mode is discussed 
second. 


Axial/Transverse Fatigue Crack Modes 

Prediction of fatigue behavior of a notched laminate utilizes 
the static failure model in conjunction with laminate property 
degradation in the notch vicinity. Predicted fatigue behavior 
includes residual strength, crack propagation, and fatigue failure 
Figure 12 illustrates the method for the axial and transverse 
fatigue crack propagation modes. 

At N = 1 cycle, all material properties are at their initial 
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static values. As cyclic loading progresses , material behavior 
at the notch will change . These changes with N can be determined 
directly from lamina properties by methods described in the pre- 
ceding subsection. The changed material properties after a given 
N are used as inputs to the static failure analysis model to 
predict the residual strength of the notched composite and the 
cprresponding failure mode. In Figure 12, static behavior at 
several decades of N computed from the static failure analysis is 
shown. In the lower left are curves of axial **crack" length a 
(inelastic length) versus applied stress a. The axial mode fail- 
ure stresses are the values of a at which the axial cracks 

extend to infinity. Overstress in material adjacent to the notch 
0 versus the applied stress o is plotted in the lower right. 

JML 

The transverse mode failure stresses are the applied stresses 

at which material adjacent to the notch fails in axial tension. 

The resulting initial static failure prediction (N = 1) in 

Figure 12 is that the notched laminate will fail by transverse 

crack propagation at a stress After ten cycles has 

occurred, a similar situation is predicted. However, the transverse 

crack propagation failure strength has now increased and the 

axial crack mode strength has decreased. The resulting residual 

( 1 ) 2 

failure stress is = 0,p . At 10 and all subsequent cycles, 

the reverse situation occurs. Either is greater than a^, as 

occurs for N = 10^, or the maximum overstress in the material 

adjacent to the notch never exceeds the unnotched laminate strength. 

2 

Therefore, for N > 10 cycles, the transverse crack propagation 
mode never occurs . 

The information contained in Figure 12 for residual strength 
and axial crack growth can be utilized as shown in Figure 13 to 
predict the residual strength, the fatigue lifetime, and axial 
fatigue crack growth with N. In Figure 13a, residual strength by 
both the transverse crack propagation mode and the axial crack 
propagation mode are plotted versus number of cycles N. As seen in 
the figure, transverse cracking will occur until approximately 20 
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cycles at which time axial crack propagation becomes the mode of 
failure.. The notched composite residual strength is always deter- 
mined by the lower of the two curves . In this hypothetical example 
the composite residual strength increases until N = 20 cycles, 
then decreases until the strength of the composite becomes equal 
to the maximum cyclic stress. The composite will fail on the 
succeeding cycle, N^, which is the resulting notched composite 
laminate lifetime. 

The static axial crack length curves of Figure 12 may be uti- 
lized to calculate the axial fatigue crack growth with number of 
cycles. Initially, the cyclic stress level is below that for which 
any axial crack will form. Somewhere between 1 and 10 cycles an 
axial crack will begin to form. This occurs when the axial 
crack initiation stress drops to the cyclic stress level. At any 
given number of cycles, the axial crack length will grow to the 
length indicated by the intersection of the maximum cyclic stress 
level with the appropriate crack length curve. The intersection 
points of the maximum cyclic stress level and crack length curves 
for different N determine the axial fatigue crack growth curve. 
Figure 13b shows the results for this example. 

Off-Axis Fatigue Crack Mode 

In principle, the same procedure used for axial crack prop- 
agation can also be applied to the off-axis fatigue crack prop- 
agation mode. Special consideration is necessary, however, to 
treat the combined stress state in the off-axis x'y' coordinate 
system. A calculation must first be made of the stress state 
, Oy, , and o^iyi around the notch root which occurs during 
fatigue loading. This is performed in a manner identical to that 
used to obtain the stress state around the notch for static loading 

With the combined stress state around the notch computed, it 
is necessary to determine how the laminate x'y* shear properties 
change under the combined cyclic loading. To do this, there are 
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two methods which may be utilized. The first method uses the 
laminate analysis technique discussed earlier to determine layer 
stresses which occur due to cf , , ^ , , and cf , , . The second 

method is completely empirical and utilizes an "effective" shear 
stress concept where the effects of combined stresses are all 
lumped into an increased value of shear stress at the notch root. 
Each method is described in more detail in Appendix E. It is 
expected that the laminate analysis technique will be the more 
accurate method. The effective shear stress method is an alter- 
native which can be used if only laminate shear fatigue data is 
available (instead of lamina date) in the xy coordinate system. 

Stresses in the Notch Region 

Determination of the stress state in the vicinity of the 
notch is critical to performing the fatigue analysis. Also, once 
the lamina stresses have been determined in the notch region, it 
is necessary to have a method for determining the effects of the 
resulting combined stress state on lamina material property 
degradation. Therefore, discussions of these topics are presented 
here. More detail can be found in Appendices, D and F, respec- 
tively. 

Laminate stress state near notch .- In a notched laminate, 
stresses and strains will vary with distance from the notch tip. 

The largest strains are expected to be at the notch, with decaying 
values as distance from the notch is increased. It would be 
beyond the resources of the present program to perform an exact 
determination of the degradation of laminate properties as a 
function of coordinates x and y. Therefore, an approximate average 
stress level in the notch root area is chosen , and it is assumed 
that all material in the notch vicinity degrades uniformly under 
this average stress state with fatigue cycling. 

For a notched laminate under a cyclic tensile stress a 

max , 

the stress state near the notch will primarily be axial tension 



and axial shear. The axial tension is due to the applied o 

max 

plus stress concentration effects, and the axial shear is due to 
the tendency of the notch core region to pull out. The tensile 
stress will exist only outside the core region of the notch, while 
the shear stress will exist inside as well. A close estimate of 
the tensile stress a” can be made by calculating the maximum 
overstress in the material adjacent to the notch, using 

the static failure model. Similarly, the shear stress t” can be 
closely approxmiated by the maximum shear stress in the notch 
region computed from the static failure model. Appendix D gives 
other, more approximate, methods for estimating a” and x”. 

Appendix D also presents an approximate analysis of the 
fatigue failure characteristics in the core region (t^ only) 
versus the adjacent overstressed material (x” plus a”). It is 
shown there that final fatigue failure predictions should not be 
significantly affected by ignoring the tensile stress o” in the 
notch region. Although stress states in some laminae can be 
different when one considers a^, the final laminate failure mode 
and failure stresses are essentially the same. Use is made of 
this result in examples given in Appendix G. 

Combined stress effects on lamina fatigue .- Interaction 
effects of stress components on static failure have been known 
and studied for some time (see, for example. References 7 and 8) . 
Experimental and analytical results indicate that the main in- 
plane stress interactions appear to be between transverse normal 
stress, ^22' axial shear stress, <^^2* reasonable to 

assume that static interaction results will be qualitatively true 
for fatigue loading. Under a laminate stress state of shear and/ 
or tension, layer stresses will generally be combined axial normal 
stress, transverse noimal stress, and axial shear stress. In the 
present formulation, the common assumption of uncoupled fatigue 
behavior between axial normal stress and the other in-plane stress 
components (such as assumed in Reference 9) will be maintained. 
Due, however, to the complexity of the experimental information 
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which is necessary to fully utilize the Reference 9 theory, a less 
complicated interaction between 022 proposed. Fatigue 

tests of [ 90 ] and [± 45 ] laminates are relatively easy to perform. 
The [ 90 ] tests will give transverse tension fatigue data directly. 
With some approximations or additional analysis work, [± 45 ] tests 
can predict the axial shear behavior of a [ 0 ] laminate. It is 
desirable to be able to predict lamina fatigue behavior under 
combined transverse tension and in-plane shear knowing fatigue 
behavior either in axial shear only ([± 45 ] laminate in tension), 
or in transverse tension only. 


Appendix F describes the development of such a theory. A 
quadratic interaction is assumed to occur between transverse ten- 
sion 0^2 axial shear a^2- physical basis for the inter- 

actions is described in the Appendix. If lamina fatigue data is 
known for either O22 or 0^2' effects of combined stresses 

can be computed from the theory. If fatigue data exists for both 
stress states, then 0^2 ” ^22 ^^^o^^otion is computed from Q22 
data when 0^2 relatively low, and from 0^2 ^^^a when 022 is 
relatively low. 


In the present work, therefore, axial tensile fatigue stress 
is assumed to cause changes in only axial Young's modulus and 
axial tensile strength. Combined effects of transverse tension 
or axial shear will be felt (through one or both of two quadratic 
interaction formulae) on axial shear modulus and axial shear 
strength, and on transverse Young's modulus and transverse tensile 
strength . 


Fatigue Analysis Procedure 

This subsection presents the methodological framework which 
has been developed to analyze the fatigue behavior of notched 
composite laminates. The analysis procedure utilizes.: the static 

failure model; lamina fatigue properties in axial tension, axial 
shear, and/or transverse tension; and constant strain laminate 


26 . 



analysis techniques. 

A notched laminate with known fiber orientation properties 
and constituent layer elastic, static failure, and fatigue prop- 
erties is under cyclic fatigue loading at a maximum applied stress 

of a and a stress ratio R. Initial static laminate properties 

max ^ 

are known either from experiment or from a laminate analysis 
prediction based on layer behavior. 

The fatigue analysis is begun by first predicting the initial 
static strength in the axial, transverse, and off-axis propagation 
modes . The lowest of the three failure stresses will be the static 
failure strength of the notched laminate. From the static analysis, 
the average axial tensile and shear stresses and at the notch 
root are computed. The laminate is analyzed with and as an 
applied stress state, and the stresses in each layer are computed. 

A small increment of cycles is chosen (say, a decade) . The 
maximum cyclic stress under axial tension, is then used with 

axial tensile layer fatigue behavior to compute the degraded axial 
tensile properties of each layer. Next, the combined stress 
effects of 022 °12 determined on the lamina axial shear 

and transverse tensile properties from lamina fatigue behavior. 

The resulting degraded properties, which will generally be dif- 
ferent for each layer, are then used in the laminate analysis to 
compute the new laminate elastic and failure properties near the 
notch. From these revised laminate properties, one can determine 
the residual strength of the notched laminate following the first 
increment in cyclic loading. 

The laminate analysis is again performed using the new layer 
properties to determine the stresses in each layer under laminate 
stresses o^ and t”. Again, the resulting layer stresses are used 
to determine the degraded properties of each layer after the second 
cyclic increment. The properties so determined are again used in 
the laminate analysis to predict laminate elastic and failure 
properties , and the new laminate properties are in turn used in 
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the static failure analysis to determine the residual strength . 

This procedure is repeated until either the residual strength 
is obtained at a given desired cycle, or until the notched lam- 
inate residual strength drops to a 

max 

This procedure is followed for the axial and transverse 
crack propagation modes, and also for the off-axis crack propaga- 
tion mode. For a given N, the residual strength will be the 
lowest of a^, a^, and The fatigue lifetime, is the number 

of cycles at which the residual strength of the laminate drops 
to the maximiam cyclic stress, 

A flow chart of this procedure is presented in Figure 14, 
Nonurtiform Fatigue Loading 

The preceding discussions of fatigue analysis have centered 
upon constant values of maximum stress stress ratio. 

The analysis can, however, be utilized to treat fatigue loading 
spectra which are not constant. It has been shown for metals that 
loading sequence can affect fatigue results, and it is anticipated 
that this effect may be even more pronounced for composites. 

Hence, the usual block loading approximation to a random spectra 
may not be a good one. However, if the loading spectrum exhibits 
tha proper behavior, it can be broken into a series of small non- 
uniform block loading sequences as an approximation to the actual 
spactrum. Provided that a nonuniform fatigue loading spectrum can 
be approximated by nonuniform block loading, fatigue of a notched 
composite can be treated using the analysis presented herein by 
tha following methodology : 

The block loading spectrum is separated into loading blocks 
bv constant a and R. Each loading block of constant fatigue 
loading parameters is then broken down into subintervals whose 
siges are dependent only upon calculation convenience. Fatigue 
laminate analysis begins on the first loading block exactly as pre- 
sented in the preceding subsection. At the end of the niimber of 


28 . 



cycles which correspond to the end of the first loading block, 

all of the current laminate and layer properties will have been 

computed. Upon beginning for analysis for the second loading 

block, the laminate now has the "degraded" properties which 

resulted from the first loading block. Provided that layer fatigue 

data is' available for the new R, the fatigue analysis can be 

continued using the new a . The analysis can be continued from 

max 

block to block until the desired residual strength properties 
have been obtained or until the notched laminate fails by fatigue 
"wearout". The key to such an analysis is obviously the avail- 
ability of layer fatigue data for many different values of R and 
and an analysis technique to estimate the effects of sequential 
loading on a unidirectional layer. 
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CONCLUDING REMARKS 


This report presents a semi-empirical analysis of fatigue of 
notched composites. The analysis incorporates heterogeneous 
material behavior and experimentally observed failure modes for 
fiber composite laminated materials. The analysis utilizes an 
approximate "materials engineering" approach along with experimental 
data . 

The overall model is composed of essentially two parts: the 

overall conceptual framework, and the details of the mathematical 
description of each physical fatigue phenomenon. Hence, the model 
is flexible since assumptions that have been made to facilitate 
analysis can be altered or replaced if desirable or necessary 
without disturbing the overall framework. The philosophy by 
which the model has been constructed is "wearout", i.e., repeated 
loading causes changes in material properties near the notch. 

These property changes result in altered notched composite residual 
strengths and eventually cause fatigue failure. 

Applications of the analysis presented in Appendix G show 
the following: 

(1) Fatigue cycling can cause an initial increase in residual 
tensile strength of notched composites. Subsequent cyc- 
ling may either cause a reversal of this trend and 
eventual fatigue failure, or it may result in a continued 
increase in strength and runout. Such behavior varies 
with notch size and laminate stacking geometry. 

(2) Static failure modes may differ from fatigue failure 
modes in the same notched laminate. Therefore, a material 
exhibiting higher initial static strength may have a 
shorter fatigue lifetime, and vice-versa. 

The fact that notched composites can have different static and 
fatigue failure modes is a result of high significance. A common 
assumption in structural reliability analysis is to assume a one- 
to-one correlation between static strength distribution and fatigue 
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strength distribution. A laminate with low static strength is 
also assumed to have a relatively short fatigue lifetime. The 
results of the present analysis indicate that in certain cases such 
an assumption may not be valid. In particular, when the mode of 
static failure differs from the mode of fatigue failure, a one-to- 
one correlation will probably not exist. The concept of static 
proof testing to ensure fatigue lifetime should therefore be 
critically re-examined. 

The static failure model has been substantiated by experi- 
ment and analysis performed in References 2 and 6. However, 
fatigue data does not currently exist which is of the form neces- 
sary to corroborate the current fatigue analysis. Therefore, a 
comprehensive program should be undertaken to verify the ap- 
plicability of the mathematical model and to correlate it with 
data from specific notched composite leuninates . A large part of 
this program would be to critically evaluate mathematical and 
physical assumptions in light of new experimental data. Since 
examples have been performed for boron/epoxy, it is logical that 
experimental correlation be attempted first with these materials. 
Testing of other fiber composites (and hybrids thereof) should 
follow provided that a first-phase testing proves the applicability 
of the model. 

The current model has distinct advantages in that only simple 
tests on simple laminates are requited as inputs to the analysis. 
These tests are: residual strength and fatigue lifetime tensile 

tests of unnotched [0] laminates, [90] laminates, and/or [±45] lam- 
inates in tension. These tests, in order to apply to different 
loading situations and different laminates, should be conducted at 
several values of and R. Although the tests required are in 

lUaX 

depth, they are uncomplicated and relatively easy to perform. 

Once the results from these tests have been obtained for a given 
material, any laminate configuration and any notch geometry can be 
analyzed. 



There are certain aspects to fatigue loading which have not 
been treated in the current model. The model is insensitive to 
stacking sequence, i.e., gross laminate properties are assumed to 
change with number of cycles and stress level, and order of layup 
does not enter the predictions. It has been observed (Reference 
2, for example) that stacking sequence can have an effect upon the 
strength and also the mode of failure during static and fatigue 
loading. In addition, it is observed for certain laminates that 
areas of debonding occur between layers as well as growth of cracks 
parallel to fiber directions. At present, it is felt that both 
of these phenomena (stacking sequence dependence and interlaminar 
failure) are caused by the same phenomenon: interlaminar shear 

behavior near the pre-existing notch or the propagating crack. it 
is therefore recommended that a program be undertaken to quantita- 
tively and qualitatively evaluate the effects of delaminations on 
fatigue behavior and provide a method for incorporating the effects 
of stacking sequence. 
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FIGURE 2. "MATERIALS ENGINEERING" MODEL FOR A NOTCHED LAMINATE 
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FIGURE 3. FIBER COMPOSITE STRESS AND COORDINATE SYSTEM CONVENTION 
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FIGURE 5. TYPICAL STATIC FA I LORE BEHAVIOR. PREDICTED FROM MODEL. 
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FIGURE 6. EQUIVALENT AXIAL CRACK PROBLEM FOR OFF-AXIS 
FAILURE 
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FIGURE 7. NOTCH ROOT UMINATE STRESS IN OFF-AXIS COORDINATE SYSTEM 



FIGURE 8. GRAPHICAL PROCEDURE FOR DETERMINING OFF- 
AXIS RESIDUAL STRENGTH OF NOTCHED LAMINATE 




Figure 9. Static Failure Stress and Mode Determination 
Procedure for Notched Composites. 
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FIGURE 10. TYPICAL FATIGUE BEHAVIOR AND LIFETIME OF 

[0] LAMINA REQUIRED TO GENERATE LAMINATE 
FATIGUE BEHAVIOR 



(a)LAMINATE AXIAL SHEAR RESIDUAL STRENGTH AND LIFETIME 



(b)LAMINATE AXIAL SHEAR RESIDUAL MODULUS 



(. 75 ) (. 50 ) 

Nf N Nf 


(ciLAMINATE AXIAL SHEAR RESIDUAL FAILURE STRAIN 

FIGURE 11. TYPICAL LAMINATE FATIGUE LIFETIME AND RESIDUAL 
PROPERTY INFORMATION - AXIAL SHEAR AT 
NOTCH ROOT 
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FIGURE 12. SCHEMATIC DIAGRAM OF METHOD FOR DETERMINING 
FATIGUE LIFETIME, RESIDUAL STRENGTH, AND AXIAL "CRACK" 
GROWTH OF NOTCHED LAMINATE. 

44 . 




(a) 


NOTCHED 

COMPOSITE 

TENSILE 

STRESS 



NUMBER OF CYCLES^ N 


(b) 


AXIAL 


"crack" 

LENGTH 


o 



FIGURE 13. TYPICAL PREDICTED RESIDUAL STRENGTH. FATIGUE 
LIFETIME. AND AXIAL "CRACK" GROWTH OF NOTCHED LAMINATE. 
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APPENDIX A 


GOVERNING EQUATIONS FOR FAILURE MODEL 


Referring to Figure 2, main text , the equations of equilibrium 
in. the various regions of shear stress transfer are as follows: 

Region of In-plane Shear Damage (Axial Crack) 

Equilibrium of the center core and adjacent overstressed 
regions in the x-direction gives; 


n k h E 

X 


d2 u. 


dx2 


0 


, (A2) 


m j h E^ 


d2 Uj 

dx2 


^ + So. 


^90^ 


m 


m 


(U^ - Uj) - 


(A2) 


where 


j = number of layers in the laminate 
jgQ = number of 90 degree layers in the laminate 
h = thickness of lamina 
E = laminate longitudinal modulus 

X 

G = axial shear modulus of the laminate 
xy 



a = gross laminate stress 
<3 

GIq = contribution of the fiber bending stiffness to the modi- 
fied shear modulus G^q of the 90 degree layer 

n = 2a = notch width. 


It is obvious that only gross properties are used and no con- 
sideration is given to the fact that the axial and shear stiffness 
actually vary along the thickness coordinate. The last term in 
equation (A2) represents the modified shear stiffness of the 90 
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degree layers in the laminate and is derived herein by assxaming 
that the fibers are clamped at two points separated by a length m 
which is the assumed region of shear transfer. The modified shear 
modulus GgQ is the sum of the shear modulus of the 90 degree 

lamina and the contribution of the fiber bending stiffness in the 
90 degree layers over a length m . As m increases, the fiber 
bending stiffness decreases and G^q approaches the conventional 

value of G^y for the 90 degree laminate. In [O/+0/9O] laminates 

with 0® i 9 i 45°, the bending stiffness of the fibers in ±6 
layers is negligible. 


Equations (Al) and (A2) are rewritten as: 

U„ 


n* j A- E 


f Lf dx^ 


= 0 


(A3) 


m* j A- E 


f Lf dx2 




m*d2 V, 


m*d2 


(Ui - U 2 ) = 0 (A4) 


where = E^/V^ , 

Vj = fiber volume fraction, 

d = 0° layer fiber spacing, 

A^ = hd; 

m* = m/d = number of fibers in overstressed intact region 
in a 0° layer 

n* = m/d = number of fibers in notch region in a 0° layer. 

The above differential equations are nondimensionalized by 
assuming 


and 


U 


0,1 


"gf 


' Vj m* \% 


u 


0,1 


d 


X 


I \f ^f 

\ "xy 



h 


C d 


where = o^/V^ . 


(A5) 


(A6) 
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Thus , we have 


n’ 


~dF' 


= 0 


(A7) 


m*d^ u. 


dC' 


- (1 + r) (Uj^ - C) = 0 


(A8) 


where r = 


^90 ®90 


k G 


xy 


The solution of equations (A7) and (A8) are, 

"O = ^1 


Ui = 5 + C 2 - e~^^) 0 < C < C 


(A9) 


1 + r ^ 

where L = ( — r — ) and C, and C~ are constants, 
m* 1 2 

The boundary conditions -gy (S = 0) =0 and u^^ (5 = 0) = 0 
have already been considered. 

Inelastic Region 


The following governing equations are derived for x-direction 
equilibrium in the inelastic zone (shear stress in shear region 
= T®, shear strain < y ) 


d2 U, 


n* j A^ E 


f Lf dx2 
d2 U, 


- 2j h T® - 2 


So So^f 


m*d2 V. '“0 “1 


(U. - U, ) = 0 (AlO) 


' 1 f) ^9 

^ ^f ^Lf + j h T® - (U, - UJ 


m*d2 Vf "^1 “O' 


Sy ^^f ^ So So^f 

m*d2 Vj m*d2 


(Uj^ - U2) = 0 


(All) 


where t® = yield (failure) stress of the laminate. 
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In addition to the proceeding nondimensionalization of the 
displacements and the coordinates , , the yield shear stress 

T° is nondimensionalized in the following manner: 


V, \ 

xy f 

'gf m* ^ 


"y • 


(A12) 


where T is the nondimensional shear stress. 


Equations (AlO) and (All) appear in nondimensionalized form as: 
d2- 


n’ 




- " 2y„ - 2r (u. - u, ) = 0 


(A13) 


d u 


m» 


1 . - 




+ T„ - r(u, - u.) - (1 + r) (u, - C) = 0 (A14) 


The solutions for equations (A13) and (A14) are : 
T. 


UQ = C- -^+C2e-^ +C^e ^ e ^ e 


"1^ 


Yi^ 


-Yi? 


Ui = C + Bf (C^e ^ e -" ) + B 2 (C^ e 


^2^ 


+ Cg e 




^ Ida 
(A15) 


where 


1,2 I 


_ 2r (m*+n*) +n* ^ 1 

^ 2m*n* “ 2m*n* 


2r(m*+n*)+n* - 8r m*n*(l+r) 


!s|J5 


ahd 1 " T;: Y 


n 2 
2r ^1,2 • 


It is evident that the above solutions are not valid for 
r = 0, as in the case of a laminate with no 90 degree layers. If 
r = 0, the solutions assume the form: 


Un = + C, C + C, 


n’ 


Ui = C + Ty + C5 e^^ + Cge^^ 


S d d a 


(A16) 
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, C^, Cg and Cg are constants of integration. 


Elastic Region 

In the elastic region, the equations of equilibrium are; 


u^ 


+ 2(1 + r) (u, - u„) = 0 


m* 


d^ Uj 


- (1 + r) (u- - u.) + (1 + r) (C - u,) = 0 


The solutions of equations A17 and A18 are; 


-e-,c -6,? 

Uq = ? + e t Cg e 


-6iC 

Uj^ = c + Bg e + Cg e 


^ a < C 


where 


:= r (t* t („.2 + n.af 

1,2 1^ m*n* m*n* J 


and 


®3,4 = 1 - n* + r) 


(A17) 

(A18) 


(A19) 


and C^ and Cg are constants of integration. In order to ensure 

finite stresses at » , terms involving positive exponentials have 
been eliminated from the displacements in equation A19. This 
boundary conditions of the problem are is nondimens ional dis- 
tance from notch to axial crack tip; and a is nondimensional 
distance from notch to beginning of elastic zone - the "inelastic 
length"); 


(C = 0) = 0 

d Uq 

(C = 0) = 0 


(A20) 

\ 
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and 



(C = 

c") 

= "0 

d "0 

dC 

(e = 

r) 

d 'in 

= <5 = 


(t = 

r) 

= (c = c+) 

d 

dC 

(e = 

r) 

d u, . 

= , ^ (5 = 

dC 

“0 

(c = 

a~) 

= Uq (C = 

dC 

(c = 

a") 

d u- . 

= ^ 0 (C = 01 ) 

dC 

“1 

(c = 

a"*) 

II 

c 

II 

Q 

d 

(c = 

a”) 

"1 + 

= TT <E - “ ' 

d Uq 

( 5-*-“ 

■) = 

d u, 

~gg'~ = 1 

Uq (5 = a 

- 

• Uj^ (C = a”) = 






> 




(A21) 


(A22) 


(A23) 


(A24) 


The boundary conditions of Equations A20 and A23 have already 
been used. Hence, there remain a total of nine boundary con- 
ditions and nine unknowns which consist of the eight constants 
Cl# C2 / C3, C4 , C5, Cg , Cy, Cg and the nondimensional yield 
shear stress ty. 
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Derivation of the Modified Shear 
Modulus for a 90 Degree Laminate 

It is argued that the shear modulus of a 90 degree laminate 
is not correctly represented by the conventional value of when 

the shear transfer zone defined by the length m in Figure 2 is"^ 
of the order 10"^ because of material heterogeneity. To show that 
this is indeed true, the actual shear modulus of a 90 degree 
laminate is derived below by considering the fiber and matrix 
phases separately. 


Consider a 90 degree lamina of width 2l and of infinite 
length which is clamped at the two edges as illustrated in 
Figure A-1. Suppose that the lamina is subjected to a differ- 
ential displacement. As a result of the deformation, the strain 
component e in the matrix and fiber will consist of the bending 
strain contribution, the shear strain as a result of this bending 
strain and the shear strain arising out of the displacement com- 
patibility at the fiber matrix interface. Referring to Figure A-2 
compatibility of displacements at the fiber-matrix interface 


require that 
h* 

^ (U 




+ - 




= c ( 




tn 


- V - ^,x) 


(A25) 


where Uy 


U 


y/X 


= displacement in the 



y direction 


Y = bending shear strain 


G = shear modulus 


The subscript f refers to the fiber and m to the matrix. It 
is assumed that Uy is independent of the y direction. Also, 
from stress equilibrium = ^xy* 

The solution for x^y from Equation (A25) is. 


^xy 


Gf 

TT+g^ 


(u 


y»x 


+ Yf 92 



(A26) 


where 


Gm 


= fiber volume fraction. 
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The total strain energy (S.E.) of the fiber and matrix consists 
of the bending strain energy, bending shear strain energy, and 
shear strain energy arising from the shear stress Tj^y . Therefore, 


S.E 


’ { '*f 


XX 


2 

) dx 


K E-2 I 2 K e_2 I 2 \ 

m / 0 


y ,xxx 


dx 


Gf h* c Gf 


kV 


(ltg2) 


^ L "y.x ^ ^2 57^' ''=■ 


(A27) 


where E = Young ' s modulus 

I = moment of inertia 
a - cross sectional area 
K = shear coefficient and 
V = shear force. 

The form of displacement assun^d is. 


( 


U 


y 


= BX' 


(A28) 


where B is a constant. The differential displacement will 
give rise to a shear force S across the cross section at 
X = i,. Equating the internal strain energy to the external, 
work done by S , 


6Ef If B^ + 18 K if- 

af Gf 


B^ i 


+ h*9B2 

gJTT+gJT 2 ® 




2 2 

+ 93 ^ + f 93 


')■ 


= s B£* 


(A29) 
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where = 


2Vf K Ef If 


The bending and bending shear strain 


energies in the matrix are neglected. 

Further manipulation of Equation (A29) yields, 


’90 


Vf Ef 


(h+2c) Aa2 2 
2 


+ 4 




* ~T~ I “I 




(A30) 


It is obvious that Ggo represents the modified shear modulus 
of the 90-degree lamina. For (h/'i)<<l, Equation (A30) can be 
written as, 


Ggo 


0.9 


(1-Vf) 



m 

5e Ji- \ 

* Of 1-Vf j 


(A31) 


which is conventional shear modulus G^y of the 90-degree 
lamina as obtained from the stiffness in series model. Note 
that the factor 0.9 results from the approximate nature of the 
displacement in Equation (A2 8). In Equation(A30)the first term on 
the right hand side may be designated as the contribution to 
the shear modulus from the bending stiffness of the fibers and 
will be termed as Ggg • Hence, 





+ IC 



(A32) 


If the shear strain t /G reaches the yield value, it is 

apparent that - Thus if the shear strain attains its 

ultimate value when fiber matrix debonding results, the 90 degree 
lamina will still have a residual stiffness G^q to resist differ- 
ential displacement. For a through-the-thickhess axial crack to 
propagate either in a [0/90 or [0/±45/90]g laminate, the fibers 

in the 90 degree lamina should break. This will be the case when 
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the maximum strain exceeds its ultimate value. Maximum 

occurs at x = I at a distance farthest from the neutral axis in 
the fiber and its magnitude is. 


e 






(A33) 


Hence the maximum fiber stress is, 

/ 



(A3 4) 


where = A£ ^ = differential displacement over the length I . 

Figure A-3 illustrates the variation of ^90 with m* , the 

number Of intact fibers in the 0° layer. It is clear that for 
m* < 5, the shear modulus underestimates the stiffness of 

the 90 degree laminate. However, for m* > 9 , approaches the 

value G 

xy. 
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FIGUr A-1, 
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FIGURE A-2, DEFORMATION OF A [90] LAMINATE SUBJECTED TO 
AN IN-PLANE SHEAR LOAD. 


so 



m 


FIGURE A-3; VARIATION OF 90° LAYER SHEAR MODULUS. 

WITH THE NUMBER OF INTACT FIBERS IN 0° LAYER, m 



APPENDIX B 


COMBINED STRESS EFFECTS ON LAMINATE STATIC AXIAL CRACKING 

Performing an exact analysis of axial cracking including 
transverse tension a„ , and applied axial shear a , , in addition 
to axial tension is difficult. Therefore, the effect of 

shear and transverse normal stresses is taken into account in an 
approximate fashion. The presence of laminate gross transverse 
tension and axial shear is treated by superimposing 

these stresses on the axial shear stress at the notch root which 
would exist due to only o^, (see Figure 7, maintained). The 
main effect of this superposition upon the shear lag analysis will 
be to increase the likelihood of x'y' axial shear failure. 

Throughout the loading process , the presence of applied axial 

shear stress o , , will reduce the amount of axial shear stress 

X y 

required from the shear-lag analysis to cause axial cracking. 

The effect of applied shear is therefore taken into account 
reducing the x'y' laminate shear strength t° ' directly by 
The effect of the transverse normal stress , is assumed to 
further reduce the shear strength according to an empirical quad- 
ratic interaction formula. 

The resulting effect of axial shear and transverse stresses 
is, therefore, to change the laminate axial shear strength in- 
putted into the standard tensile strength analysis. This changed, 
or effective, shear strength t®' is found from the equation 
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where 


to . 

e 


= the effective shear strength for given o 


y ' 


and 0 


x'y’ 


a 


P 

y' 



+ o 


t 

yl 


2 


o 


d 

y' 




tensile strength of the unnotched laminate in the 
y’ direction 



compression strength of the unnotched laminate in 
the yV direction 


• d _ 


• O I _ 


x'y ' 




T®' = shear strength of the unnotched laminate in the 
x'y' coordinate system. 


For a given laminate and fiber orientation, increased applied 
stress a will cause proportional increases in the stresses , and 
a^,y,. This means that the effective, reduced shear strength, t ® * 
will change (and generally decrease) with increased tensile load- 
ing. 
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APPENDIX C 


STATIC FAILURE ANLAYSIS EX7VMPLES 

This appendix gives examples of the notched laminate static 
failure analysis procedure- Unidirectional and multidirectional 
laminates are given separate consideration to bring out details 
of the analysis and to illustrate the major differences between 
the two systems- Laminates chosen for example computations are 
[ 0 ] and [ 0 2 / 1 4 5 ] boron/epoxy - 

Unidirectional Systems 

Unidirectional composites are special from the failure sense 

in that the only crack propagation modes are along the zero degree 

fiber in the load direction or transversely across the composite 

normal to the loading direction. As such, only the axial and 

transverse crack propagation modes need to be considered, and off- 

axis cracking is not a possibility. Input data to the failure 

analysis are the unnotched laminate unidirectional tensile strength 

the unnotched laminate shear strength t°, the laminate shear 

modulus G , the laminate extensional modulus E , the laminate 
xy ' x' 

ultimate strain in axial shear T^2.t' size 2a, and over- 

stressed zone size m. 

As an example of the procedure for determining static failure 
of unidirectional composites, a boron/epoxy [0] laminate is analy- 
zed. As mentioned previously, the analysis predicts the stress 
at which an axial crack grows to infinity under an applied ten- 
sile stress, O- The analysis also determines the maximum stress, 
which! exists in the material adjacent to the notch and the 

SCFM 

corresponding applied laminate stress, a^, at which it occurs. In 
all cases, 

It is to be expected that failure strengths of a notched com- 
posite will vary with the constituent material properties of the 
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laminate and the notch size. Figure C-1 illustrates, for a typical 
[ 0 ] boron/epoxy composite laminate with a 2.54cm notch, how the 
axial failure stress and the adjacent material overstress with 
corresponding laminate tensile stress can vary with properties. 

In this example, laminate shear failure stress, , and the shear 
failure strain, Yult' varied and the shear modulus is held 

constant. Generally speaking, lower values of x° and 
give notched laminates with lower axial crack propagation failure 
strengths. Also, this same trend of lower t® and reduce 

the laminate stress concentration and therefore lower the tendency 
for a transverse crack to develop. For increasing values of x® 
and Yult' reverse situation is true. Axial crack propagation 

in the axial direction is impeded with increasing x® and Yuit' 
and stress concentrations in material adjacent to the notch are 
increasing. Therefore, tendency of the notched laminate to fail 
by a transverse crack propagation mode is enhanced by increasing 
x“ and Yult- 

Material properties for a typical [0] unnotched B/Ep laminate 

are given in Table C-1. To obtain predictions of the failure model 

for static notched laminate strength, the typical failure shear 

2 

stress and shear failure strain from Table C-1 of 105 MN/m and 

.038 m/m, respectively, are entered into Figure C-la. The resulting 

prediction for axial crack propagation failure stress is about 380 
2 

MN/m . In Figure C-lb, the predicted value of maximum composite 

2 

overstress is cTscfm~^^*^ MN/m , which is considerably less than the 

2 

unnotched tensile strength of 1300 MN/m for the laminate. The 

transverse tensile failure mode will, therefore, not occur. Static 

failure predictions for this laminate and notch size are that the 

notched composite will fail due to axial crack propagation at a 

2 

stress level of 380 MN/m . 


Figure C-2 illustrates the variation of axial crack propagation 
mode failure stress and maximum composite overstress for a 0.635 


The model predicts = 580 MN/m^ 


cm notched boron/epoxy material . 

2 

and ~ 1030 MN/m . Since the composite overstress is still 
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2 

less than 1300 MN/m (the failure stress for the unnotched lamin- 
ate) , transverse crack propagation will still not occur for this 
laminate. However, the 0.635cm notch axial crack propagation mode 
strength is considerably higher than for the 2.54cm notch. This 
trend generally follows observed experimental behavior for the 
effect of notch size on strength of notched composite laminates 
(Reference 2) . 


Multidirectional Systems 

Failure predictions for multidirectional composite laminates 
differ significantly from [0] composite laminates in two major 
ways. First, as discussed previously, there are three possible 
failure modes rather than two. These are; axial crack propaga- 
tion along a zero degree fiber, transverse crack propagation a- 
cross the specimen normal to the direction of loading, and off- 
axis crack propagation along one of the off-axis fiber directions. 
The second main difference is in the stress-strain behavior in 
axial shear. For [0] laminates, the stress-strain behavior in 
axial shear is nonlinear, and can be approximated by a bilinear 
elastic/plastic stress strain curve. For multidirectional lam- 
inates , the constraining effect of the fibers will tend to elimin- 
ate matrix nonlinear effects. As a result, stress-strain behavior 
for almost any loading condition on a multidirectional laminate 
will be nearly linearly elastic to failure. Therefore, to a good 
degree of approximation, multidirectional 

laminates, therefore, shear stress-strain behavior can be completely 
described knowing only t° and tO] laminates, however, 

Yuit^ T° , and can all be independent variables. 

As an example of the method used to predict the failure mode 
and failure stress of a notched multidirectional laminate, a boron/ 
epoxy [ 02 /^ 45 ] laminate with a 2.54cm slit center notch is selec- 
ted. It is assumed that this laminate has the overall properties 
shown in Table C-2 which are predicted from unidirectional proper- 
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ties given in Table C-1. The laminate predictions were made using 
a s.tandard constant-strain laminate analysis technique and uti- 
lized a maximum stress failure criterion. While agreement with 
experimental data is not exact, the predicted values to be used 
are reasonably close to experimental values. 

Axial/transverse crack propagation .- Laminate properties 
shown in Table C-2 were utilized in the analysis described in 
Appendix A and the main text. The variation of and '^ 3 Qp]\,j with 
the laminate shear properties t® and were calculated from the 
analysis. Figure C-3 shows the results in graph form. The fail- 
ure stresses in the axial and transverse crack propagation modes 
can be determined from Figure C-3 if laminate shear properties are 
known. It is seen that the predicted value of axial crack prop- 
agation stress for the specific laminate properties in Table C-2 
2 

is 0- = 250 MN/m . From Figure 11b, the maximum overstress in 

” 2 
material adjacent to the notch is <^gQpj^ ~ MN/m . When this 

number is compared to the 710 MN/m^ tensile failure stress of the 

unnotched laminate, it is concluded that transverse crack propag- 

agation will not take place. This notched laminate is therefore 

2 

predicted to fail by axial crack propagation at = 250 MN/m . 

Off-axis crack failure .- Once the failure stress has been 
determined for the axial/transverse modes , it is necessary to de- 
termine the stress which would cause off-axis cracking. In this 
case, off-axis cracking could occur along fibers in the ±45 degree 
directions. For the off-axis analysis, the laminate is considered 
in an x'y' coordinate system where x' lies along the -45° fibers 
in the direction of the crack. The resulting laminate is a 
[0/+4 52 / 90 ] laminate in the x'y' coordinate system as shown in 
Figure C-4a. The laminate stress state in this coordinate system 
is = a/2, = a/2, and = °/2 • 

For computational purposes, it is assumed that a [0/+4 52/90] 
laminate exhibits the same notched composite failure characteris- 
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tics as does a [0/±45/90] laminate. Elastic and unnotched failure 

properties were predicted from a constant strain laminate analysis 

using a maximum stress failure criterion. Results for these are 

shown in Table C-3 for [02/^45], [Oy-^^S^/^O] , and [0/±45/90] boron/ 

epoxy laminates. For comparison, experimental results for the 

[ 02 /^ 45 ] laminate from Reference 10 are also shown in Table C-2. 

Note that only the axial shear modulus, G , and the shear failure 

xy 

stress, T°, differ to any great extent between the 

and [0/i45/90] laminates. Therefore, results of the failure model 
for a [0/±45/90] laminate are used with the shear modulus and 
failure shear stress for the [0/+452/90] laminate to predict the 
X* axial failure stress, It is noted that since we are deal- 

ing with a slit notch, the effective notch width in the x*y' co- 
ordinate system is now 2.54cm x cos 45®, or 1.80cm. (Figure C-4a). 

The off-axis propagation mode strength for the [O 2 /- 45 ] lam- 
inate is given as a function of and * in Figure C-4b. 

Using the shear modulus and axial shear strength for the [0/+45<~/90] 

2 ^ 

laminate gives a value of 445 MN/m for o \ (not shown in the 

A 

figure) . This value, however, is not the failure stress of the 

notched laminate in the x' coordinate system, because it includes 

the effects of only o , = a/2. The effects of a , and a , , have 

X Y X y 

yet to be determined. As evidenced by equation (Bl) , Appendix B 
the failure stress in x'y' shear to be used in the analysis is 
dependent upon the values of a^, and this laminate, 

0 , = 0 , = a , , = 0 , . Equation (Bl) may therefore be rewritten 

X y X y 

for this particular system to read: 


. o • _ 


^A' 


Oa. t 20 ^, 0 ^, 


= 1 


(Cl) 
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This equation may be solved for t'” if a , is known. However, the 

particular value of desired is the one which corresponds to the 

failure stress a which is at present unknown. Therefore, following 

the procedure outlined in the main text, several values are chosen 

2 

for between 0 and 445 MN/m , and the corresponding values of 

are calculated from equation (C-1) . The curve of Figure C-4b 

is then entered with each of these values and a corresponding 

is computed. The calculated which coincides with the 

value of the assumed a , is the desired failure stress. In this 

^ 2 

case, the failure stress is found to be 176 MN/m . The correspond- 
ing a®* is shown in Figure C-4b. Since for 0 = 45° cFtv i = cr i = cr/2, 

the axial tensile stress in the x-direction which causes 45® off- 

2 

axis cracking is a = 352 MN/m . 

Predicted failure mode and failure stress .- The stress at 

which off-axis cracking occurs is higher than the failure stress 
2 

of 250 MN/m predicted for the 0° axial crack propagation mode. 
Therefore, the static failure mode predicted for the 2.54cm notched 

B/Ep [0-/±45] laminate is axial crack propagation in the x-direction 

^ 2 
at an applied tensile stress of 250 MN/m . 
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TABLE C-1 


TYPICAL [0] BORON/EPOXY ELASTIC 
AND FAILURE PROPERTIES 


MODULI (GN/m2 ) 

INITIAL 

TANGENT 

SECANT TO 
FAILURE 

Young's 

207 

207 

Young's £^2 

22.8 

17. 4 

Shear Modulus Gj ^2 

6.4 

2.8 

Axial Poisson's 
Ratio, v-^2 

0.21 

0.21 


FAILURE STRESSES (MN/m^) 

Axial Tension 

1296 

Axial Compression 

2496 

Transverse Tens ion o 2 2 ^ 

69 

Transverse Compression ^ 22 ^ 

310 

Axial Shear t °2 

105 


FAILURE STRAINS (m/m) 

Axial Tension 

6.27 X 10“3 

Axial Compression 

12.07 X 10”3 

Transverse Tension 

3.95 X 10“3 

Transverse Compression 

17.80 X 10-3 

Axial Shear (Engineering) 

38.25 X 10-3 






BORON/EPOXY LAMINATE ELASTIC AND FAILURE PROPERTIES 
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* Constant strain analysis , max. stress failure criterion 




LAMINATE 

AXIAL 

SHEAR 

STRENGTH 

T° 

MN/m2 



ULTIMATE AXIAL SHEAR STRAIN. Yu|| . M/m 

(a) NOTCHED LAMINATE TENSILE STRENGTH. a 
AXIAL CRACK FAILURE MODE. 



ULTIMATE AXIAL SHEAR STRAIN, m/M 

( b) COMPOSITE OVERSTRESS IN MATERIAL ADJACENT TO 
NOTCH, AND CORRESPONDING FAR FIELD 

APPLIED GROSS STRESS, 


FIGURE C-1. STATIC FAILURE BEHAVIOR OF [0] BORON/EPOXY NOTCHED 
LAMINATE. NOTCH SIZE = 2.5^1 cm. 
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ULTIMATE AXIAL STRAIN, Yuit-' M^M 

(a) NOTCHED UMINATE TENSILE STRENGTH, o, , 
AXIAL CRACK FAILURE MODE ^ 



(b) COMPOSITE OVERSTRESS NEAR NOTCH AND 

CORRESPONDING FAR FIELD APPLIED GROSS STRESS, o^, 

FIGURE C-2. STATIC FAILURE BEHAVIOR OF 101 BORON/EPOXY 
NOTCHED LAMINATE. NOTCH SIZE = 0.635 CM:. 
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FIGl'Rl! C-3. STATIC FAILURE BEHAVIOR OF BORON/EPOXY [02/±A51 NOTCHED LAMINATE. 

NOTCH SIZE = 2.5A CM. 
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SLIT NOTCHED LAMINATE. NOTCH SIZE = 2.54 cm. 


APPENDIX D 


NOTCH ROOT STRESS STATE CONSIDERATIONS 

Knowing the notch room stress state during cyclic loading is 
critical to determination of material property degradation in the 
notch region. This appendix discusses methods for obtaining these 
stresses, evaluates the relative importance of laminate axial shear 
and axial tension, and discusses the relative effects of lamina 
stresses on laminate property degradation. 

Laminate Stress State Near Notch 

Inside the core region near the notch, laminate material will 
be under a stress state which is primarily axial shear. Outside 
the core region in the overstressed material adjacent to the notch, 
axial tension will be present in addition to axial shear. In 
general, the stress levels will decrease with increasing distance 
from the notch. The static failure model is capable of predicting 
the variation of shear stress, t, and the composite overstress in 
adjacent material, as functions of distance x from the notch. 

(t is assumed constant with y in the shear transfer region, and 

is assumed constant with y in the overstressed zone) . The- 
oretically, these stresses could be used to predict material prop- 
erty degradation as a function of x and number of cycles N. A 
static failure analysis could be developed which allowed material 
properties to vary with x, and a fatigue analysis could be per- 
formed taking into account the spacial variation in fatigue- 
degraded properties . 

Such a procedure is beyond the scope of the present effort. 
Instead, it is assumed that enough accuracy can be obtained by 
selecting an average stress state in the notch region, and allowing 
all material to degrade uniformly. This average stress state 
will be one of axial shear, t^, and/or axial tension, 
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The most accurate method of predicting and under an 
applied cyclic stress is to utilize the static failure an- 

alysis under the laminate stress o . The maximum shear stress, 

ITld.X 

and maximum overstress in adjacent material, t>e 

computed from this analysis. It can therefore be assumed that, 
the notch root stresses causing the major material property 
changes are: 

= ^1 (Dl) 

= ^2 '^SCFM 

Where the k's are empirical constants which can be chosen either 
analytically to represent "average" stresses or experimentally to 
fit test data. 

Stress concentrations predicted by the static failure analysis 
are typically 1.2 - 1.5. A second less accurate method of select- 
ing notch root tensile stress would be to set: 

= k 3 , k 3 = 1.2 - 1.5 (D2a) 

A companion assumption for estimating average notch root 
shear stress is to assume that the maximum shear stress at the 
notch bears the same relationship to laminate axial shear failure 
strength as the laminate cyclic maximum stress, ^j^ax' 
to the initial static axial mode failure strength of the notched 
laminate, * Since this is approximate, it is anticipated that 
an empirical factor might be necessary to achieve the necessary 
correlation with experimental data. This results in the following 
equation: 

•T = k^ , . max ,^0 (D2b) 

A 

An alternate estimate of to go with (D2a) is the reduced lam- 
inate axial shear strength, x® ^ which would cause static fail- 
ure of the notched laminate to occur at 

= k T° 

^5^amax (D2c) 

in all cases, the k^^ are empirical constants of order unity. 
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Effects of Stress Components in Laminate Fatigue 

It is desirable to determine the relative importance of lam- 
inate axial shear stress, and axial normal stress, a^, to the 
fatigue degradation of material around the notch. If does not 
contribute significantly to degradation, then the core and over- 
stressed regions will degrade to approximately the same degree. 
Fatigue analysis can be performed using a simple shear stress 
state, which is easier to handle analytically than combined tension 
and shear. Also, the lamina fatigue data base generally required 
for analysis includes fatigue behavior under axial normal stress, 
transverse normal stress, and axial shear stress. If it could be 
demonstrated that one of the stress components example) 

caused negligible damage compared with the others, the required 
data base for fatigue analysis could be reduced. 

To investigate the above topics, selected boron/ epoxy laminates 

from the family [0 . /-45 ./90, ] are analyzed. Standard constant 

1 3 

strain laminate analysis techniques are used to determine layer 
stresses under laminate stress states which approximate the states 
of stress near the notch. Two stress states are considered: (1) 

unit xy shear stressing and (2) equal unit tensile stressing in 
the X direction coupled with unit xy shear stressing (see Figure 
3 main text) . 

Laminates analyzed are [0] , [0/90] , [±45] , [0/±45] , [0 2/±45] 
and [0/±45/90] . The results are shown in Tables D-1 andD-2. The 
stresses O 22 ^12' column 4, Table D-2, are the 

layer stresses due to the unit laminate gross state stress shown 
in the far left hand coliimn. is normal stress in the fiber 

direction, O 22 normal stress transverse to the fiber direction, 
and aj^2 axial shear stress. Figure 3b). The relative sizes of 
these stresses are important to deciding whether or not a partic- 
ular stress component causes damage. Accordingly, their ratios 
are shown in columns 5 through 7. A lower case t indicates 
normal tensile stress and a lower case c indicates normal compres- 
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sive stress. Since the sign of the shear stress is irrelevant to 
failure, it is always positive. The symbols t/c, etc. in column 
6 indicate that is tensile and O 22 compressive, etc. 

As a rough gauge of the importance of a particular stress 
component to the fatigue degradation occurring in a particular 
layer, the relative stress ratios for each layer are compared to 
ratios of the static failure stresses for a unidirectional layer. 
The static failure stresses for boron/epoxy are shown in Table 
C-1, Appendix C, and their ratios are shown in Table D-1. The 
computed stress ratios for each layer and loading condition are 
shown in columns 5-7, Table D-2. In column 8, Table D-2, is 
given an estimate of the fatigue failure mode for each layer 


which was obtained by comparing each layer stress ratio with the 
corresponding failure stress ratios. As an example of the proced- 
ure followed, consider the f±45] laminate under the pure shear 
stress state a =1. Laminate analysis shows that axial tensile 
and transverse tensile stresses were both enormously higher than 
axial shear stresses. Therefore, axial shear is not a likely 
fatigue failure mechanism with the [±45] laminate. In the +45® 


layer, the ratio of to 022 is 14.6, with tensile and 
compressive. The ratio of equivalent static failure loads. 


^22 


shown under t/c as 4, is obviously much lower than the relative 


stress ratio in the laminate. This would indicate that axial 


tensile stress is the most likely candidate for causing fatigue 
failure in the +45® layer. On the other hand, the -45® layer 
experiences axial compression and transverse tension. The stress 
ratio of 14.6 c/t, when compared with a static failure ratio of 
36 c/t, indicates that the tensile stress is probably not high 
enough to cause significant fatigue damage, and the main fatigue 
damage most likely results from a matrix damage mode in transverse 


tension. 


Based upon the individual layer results, an estimate is made 
of the layers in each laminate which are most likely to fail under 
xy fatigue loading. A prediction is then made of the mode of fail- 
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ute. Results are shown in column 9 of Table D-2. 


A camparison of results for all laminates analyzed under the 
pure shear and combined tension shear stress state shows the 
following: 

(a) Significant differences in lamina stress states exist 
between pure shear and combined tension/shear cracking, 
but 

(b) These differences do not appear to significantly affect 
the outcome of fatigue failure predictions except pos- 
sibly for the [±45] laminate. 

It is concluded that inclusion of axial laminate tension in the 
notch root stress state for fatigue analysis may not be necessary. 

A final judgement, however, should await more complete analysis 
and experimental data. 

It is seen from column 9, Table D-2, that all lamina stress 
components play an important role in fatigue behavior in the notch 
region. Lamina fatigue data base tests should therefore be per- 
formed under axial tension, transverse tension, and axial shear. 

It is noted that axial shear is important only to unidirectional 
[0] and bidirectional [0/90], [±45] laminates. For typical aero- 

structure layups such as [0/±45] , [0/t45/90] , and [02/^45] , shear 
is not an important source of lamina degration. This indicates 
that, for these laminates, lamina fatigue data base tests need 
only be performed on tensile specimens with [0] and [90] configura- 
tions to develop and C 22 
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TABLE. D-2. PREDICTED IMPORTANCE OF STRESS COMPONENTS TO FATIGUE OF NOTCHED LAMINATES 
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tension, TT = transverse tension, S = axial shear, 
shown are stress component ratios between inodes most likely to fail. 



APPENDIX E 


OFF-AXIS LAMINATE COMBINED STRESS STATE IN FATIGUE 


In an off-axis fatigue failure analysis, transverse laminate 
stress a , and axial shear stresses o , , are present in addition 

y y 

to axial tension a^, . This appendix presents two methods for 
determining the combined stress effects on the notch root material 
fatigue behavior. One is a laminate analysis technique, and the 
other is an empirical method. 

The laminate analysis technique of determining off-axis x'y' 
laminate shear property degradation utilizes the laminate stress 
state around the notch (Figure 7, main text) directly in a laminate 
analysis to obtain individual layer stresses. These stresses are 
used with individual layer fatigue data to predict changes in indi- 
vidual layer properties. Laminate analysis is then used to re- 
construct laminate shear properties from the changed layer prop- 
erties. It is expected that this first method will be the more 
accurate . 


The empirical method is based on the assumption that laminate 
transverse tension and axial shear can be treated by adding to the 
shear stress which would exist due to only axial a^, tension. 

• • • n ^ 

First, it IS assumed that the tensile stress a and the shear 
stress which would be present at the notch root due to the tensile 
stress can be chosen by methods previously presented in Appendix 

D. The applied’ laminate shear stress, a , , , is assumed to add 


n' 


x'y' 


Transverse tension , is assumed to cause an 


directly to x 

artificial increase in shear stress through an empirical quadratic 
interaction. The result is an off-axis effective notch boot shear 
Og which empirically takes into account the combined sttess 
state at the notch. This effective shear can be found from the 
equation : 
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la, 0^ , ^ 

2 

/ n' 

y' + y' 1 


f e 

1 "S. 1 


1 I 


= 1 


(El) 


Where t 


n’ 


effective laminate notch shear stress, and 

I 


.P 


[1 - (0y,/0P,)']^ 


0°, - a^. 

y 




+ a 


t 

y: 


2 


t c 

IS the laminate tensile strength in the y* direction, 

the compressive strength of the laminate in the y* direction, 
n * 

Once T has been determined, an x* axial fatigue analysis is 
formed using methods previously outlined. 


is 

per- 
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APPENDIX F 


INTERACTION EFFECTS OF LAYER STRESS COMPONENTS IN FATIUGE 

Interaction effects of stress components on static failure 
have been known and studied for some time (see, for example. Ref- 
erences 7 and 8) . Results of research in this area have been 
along two lines. One is empirical wherein forms for failure inter- 
action equations under combined stresses are assumed and fit to 
test data (Reference 7, for example) . The second consists of micro- 
mechanical anlaysis approaches which take account of known hetero- 
geneities and failure modes to analytically arrive at failure inter- 
action behavior (Reference 8) . Experimental and analytical results 
indicate that the main in-plane stress interactions appear to be 
between transverse normal stress, <^22' axial shear stress, 0 ^ 2 ' 

It is reasonable to assume that static interaction results 
will be qualitatively true for fatigue loading. Hashin and Rotem 
(Reference 9) have developed a quadratic interaction theory which 
empirically treats the fatigue of a unidirectional layer under 
combined axial shear and transverse normal stresses. In their 
theory, axial stress is assumed to be uncoupled from the other two 
stress components in fatigue. In the present formulation, the 
assumption of uncoupled fatigue behavior between axial normal stress 
and the other in-plane stress components will be maintained. Due, 
however, to the complexity of the experimental information which is 
necessary to fully utilize the Hashin-Rotem theory, a less com- 
plicated interaction between and a ^2 will be proposed here. 
Fatigue testa of [90] and [±45] laminates are relatively easy to 
perform, and, with some approximations or additional analysis work, 
will reasonably predict the transverse tensile and axial shear 
behavior of a [0] laminate. It is, therefore, desirable to be able 
to predict lamina fatigue behavior under combined transverse tension 
and in-plane shear knowing fatigue behavior either in axial shear 
only ([±45] laminate in tension), or in transverse tension only. 
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The physical basis for the proposed empirical interaction is 
as follows. Failure under axial shear and transverse normal stress- 
es are both related to matrix or fiber-matrix interface failure. 

If one makes the further assumption that failure under both stress 
components is phenomenologically similar on an atomistic level, 
then one could in principle predict fatigue failure properties in 
axial shear from those in transverse normal stress and vice versa. 

In the succeeding development, it will be assumed that axial 
shear 0^2 the stress component for which experimental results 
are available, and the effects of 022 upon this information is 
desired. The development for the effects of known 022 

behavior is identical, and is not presented here. 

Suppose that the complete fatigue lifetime and residual 
strength behavior of a unidirectional lamina in axial shear is 
known from experiment. From this information, it must be possible 
to determine the fatigue lifetime, residual shear strength, resid- 
ual shear failure strains at a given N for a given maximum cyclic 
shear stress level. The effect of the additional transverse normal 
stress component will be assumed to accelerate fatigue degradation. 
This will be done through an increased "effective" shear stress 
level and a resulting "effective" S^2 ^ ^ 12^^12 (here, ^£2 
lamina static axial shear strength ). It is necessary that ^^2' 
the effective S^2 he equal to the actual S^2 ^hen 022 is zero. 

Also, as 022 approaches its failure value in either tension or 
compression, it is reasonable to assume that its effect is similar 
to cycling at an effective shear stress close to the static failure 
shear stress, x°. 

Two possible forms of interaction between ^ud 022 

shown in Figure F- 1 . One form of interaction is linear and the 
other is quadratic. The quadratic curve generally gives more 
weight to higher levels of transverse stress and less weight to 

0 

the lower levels of transverse stress. In addition, it lowers 
when small to moderate amounts of transverse compression are ap- 
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plied. These are physically logical possibilities because moderate 
compressive stresses can retard crack propagation in either Mode 1 
or Mode 2 . The difficulty with this quadratic relationship is that 
for simple forms of a quadratic curve, it can predict a negative 
value. This situation is illustrated in the figure. 


The approach taken here is to utilize a quadratic interaction 
formula with a zero minimum for Sj^2* equation to be used is 

as follows: 



(FI) 


(If 


’12 


< 0 from above, set S®2 = P) 


where 


22 


22 


+ a 


22 


_c t 

°22 " ^^22 


22 


;P 


"l2 - 1 




C t • • 

and the and refer to static failure stresses in compression 
and tension, respectively. The calculated value of S®2 is used 
in place of S^2 predict layer axial shear property degradation 
from experimentally determined layer fatigue behavior. It is 
expected that the combined stress state will also cause degradation 
in transverse tensile properties. Specifically, transverse 
normal strength and transverse Young's modulus are anticipated 
to exhibit the greatest changes. It will be assumed that the 
degradations which occur in transverse normal strengths and trans- 
verse Young's modulus occur to the same degree as axial shear 
strength and axial shear modulus. Specifically, the same fractional 


reduction will be assumed in E22 G 


’ 12 ' 


and in 022 x 


12 ' 
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The procedure for determining the effect of transverse normal 
property degradations in a unidirectional layer is, therefore, to 
(a) determine, from laminate analysis techniques, the values of 

^22 ^ given layer, (b) compute the effective non-dimen- 
sional maximum stress S^2 equations (Fl) , (c) utilize an 

experimentally determined axial shear fatigue wearout curve to 
determine the change in axial shear failure strength and axial 
shear modulus with number of cycles N, and (d) calculate the de- 
gradation of transverse Young’s modulus and transverse normal 
strength by assuming that the reductions in these values are frac- 
tionally the same as have occurred in their shear mode counterparts. 

If the situation is reversed, ie, transverse normal stress 
fatigue results are available and it is desired to determine the 
effect of exactly analogous procedure- and equations are 

used. 

It is expected that the predicted effect of 022 ^22 

be most accurate when 022 is relatively small, and that the effect 
of 0^2 ^22 most accurate when 0^2 is relatively small. 

In either situation, it is expected that prediction of pure trans- 
verse tension fatigue behavior from pure shear test results and 
vice versa may not give accurate results. It is therefore recom- 
mended that both transverse normal stress and axial shear data be 
obtained if it is known that both stress components will be impor- 
tant to fatigue degradations. Equation (Fl) can then be used for 
situations where 0^2 dominant, and its counterpart for 0^2 (not 
developed here) can be used when 022 is dominant. It is noted 
that results of Appendix D indicate that 022 is dominant for common 
[ 0 ./t 4 5 ./ 90 , ] boron/epoxy laminates . 

1 J K 
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FIGURE F-1. "EFFECTIVE" SHEAR STRESS IN FATIGUE DUE 
TO TRANSVERSE TENSION. 


APPENDIX G 


FATIGUE FAILURE ANALYSIS EXAMPLES 


Since fatigue analysis of unidirectional and multidirectional 
composite laminates differs in several respects, they will be 
treated separately through examples in more detail in this appendix. 

Unidirectional Systems 

As with static failure, the fatigue failure behavior of uni- 
directional laminates differs from multidirectional laminates. 

The fatigue failure modes for a unidirectional laminate can only 
be axial crack propagation or transverse crack propagation. Off- 
axis modes do not exist because of the lack of off-axis fibers 
along which the crack might propagate. In addition, there is 
currently no reason to postulate a relationship between x°, 
and done for the nearly elastic multidirectional 

systems. Most importantly, the fatigue behavior of unidirectional 
systems will be dominated by matrix behavior since the laminate 
has no fiber systems running across an axial crack. Therefore, 
the fatigue analysis can be performed utilizing only fatigue be- 
havior in axial shear of the unidirectional layer. 

Example . - The procedure to be followed in analyzing unidirec- 
tional systems is exemplified by treating two hypothetical fatigue 
property variations of boron/epoxy laminate with a 0.64cm notch 
under cyclic loading. Figure G-1 illustrates these variations, 
termed Material A and B, respectively. Material A consists of a 

constant shear strength t°, a constant elastic shear modulus G , 

xy 

and a decreasing ultimate shear strain number of cycles 

N. Material B, on the other hand, has constant G , decreasing 

xy 

T° and increasing with number of cycles. Figure G-2 illus- 

trates how the axial and* transverse failure stresses vary with 

changes in laminate axial shear properties. Since G is constant 

xy 

for both materials, the only variables are t® and 
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Consider, first, the axial crack propagation mode, Fi^re G-2a. 
Material A and material B property variations with nuit±>er of cycles 
are shown as straight lines with arrows indicating the direction of 
increasing N. Material A exhibits a decrease in the axial fatigue 
strength of the notched composite, while material B exhibits an 
increase in the axial crack propagation strength. Assume that in 

2 

both cases the maximum cyclic stress for the laminate is 400 MN/m . 

3 

From the figure it is seen that at somewhat less than 10 cycles, 
material A strength has been reduced to the maximum cyclic stress, 
and the material fails. Material B on the other hand, exhibits an 
always increasing axial fatigue propagation strength, so that run- 
out is predicted in this case. 

It is necessary to also consider the transverse crack propa- 
gation mode. Figure C-2b shows variations of the overstress adja- 
cent to the notch the corresponding laminate tensile 

stress with laminate material shear properties. For material A, 

it is seen that the overstress is initially near L,000 MN/m^ while 

2 

the notched tensile strength is almost 1,300 MN/m . As cycling 
progresses, the maximum overstress in adjacent material decreases, 
removing any possibility that material A will fail in a transverse 
crack propagation mode. 

For material B, the initial static strength prediction is, of 
course, identical to material A: transverse crack propagation does 

not occur initially because the overstress is less than the un- 
notched laminate strength. As cycling continues, o for mate- 

rial B decreases with N, although not quite as rapidly as for 
material A. For material B, therefore, the transverse crack propa- 
gation mode never occurs, and all fatigue and residual failures 
will be due to axial crack propagation. The resulting residual 
strength versus number of cycles for material B is seen to increase 
as shown in Figure C-3b. 

The computation of axial "crack" length a, with number of 
cycles is illustrated in Figure G-4 for both materials A and B. 
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The curves of axial crack length versus applied stress are plotted 

2 3 

for properties at 1, 10 , and 10 fatigue cycles. As discussed in 

the main text, the axial crack length at any given number of cycles 

is obtained from the intersection of the ot versus o curve with the 

line representing the maximum cyclic stress. The resulting axial 

crack length growth with number of cycles is shown in Figure G-5. 

2 

Material A exhibits very little crack growth until after 10 cycles. 

It then develops a fast running crack as the number of cycles 
3 

approaches 10 . At this point, the axial crack propagates to 
infinity, causing failure. For material B, stable crack growth 

3 

is observed out past 10 cycles. 

To facilitate comparison of effect of notch size in [0] B/Ep 
laminates, static failure analysis results were plotted for a 2.54 
cm notch. The results are shown in Figure G-6, A direct comparison 
point by point for given sets of laminate shear properties show 
that larger notch size generally weakens the strength of the mate- 
rial in the axial crack propagation mode and makes failure by 
transverse crack propagation considerably less likely. These 
results are true both for static failure and fatigue failure. 

Discussion .- It is emphasized that the preceding examples of 
unidirectional fatigue behavior are examples only, and the material 
behavior in axial shear which has been utilized in the examples is 
hypothetical. It must not be taken as even an estimate of the 
actual laminate behavior. The example does serve to illustrate 
that different fatigue behavior for notched composites is predictable 
from the mathematical models. It is also noted that the variations 
of a^, and with axial shear material properties (Figures 

G-2 and G-6) are valid predictions. They provide a method of com- 
parison of notched composite behavior for different notch sizes, 
as discussed previously. 

Multidirectional Systems 

The primary region of material degradation in a notched com- 
posite laminate is the region immediately surrounding the notch. 
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In the example presented here, the axial shear stress is assumed 
to dominate the notched composite behavior, and the axial tensile 
stress is neglected (Appendix D) . Following a short discussion of 
the differences between unidirectional and multidirectional laminate 
behavior, an example is then given for fatigue behavior of a notched 
boron/epoxy [02/±45] laminate using assumed unidirectional layer 
material behavior . 

Differences between unidirectional and multidirectional 
laminates . - In addition to failure mode differences between [0] and 
multidirectional laminates there are major differences between 
the two types of laminate which are important from a fatigue stand- 
point. These differences relate mainly to the kinds of layer fati- 
gue and failure information necessary to predict laminate fatigue 
behavior. 

A [0] laminate under axial tension has no fibers extending 
across the axial crack which is growing from the notch. The be- 
havior of this laminate in axial shear, therefore, is primarily 
matrix controlled (although the fibers contribute by their presence 
through geometrical factors) . Since unnotched 0® specimens gener- 
ally exhibit excellent fatigue characteristics, degradation of a 
unidirectional notched laminate due to the tensile stresses present 
will be minimal, and axial shear fatigue behavior will dominate. 
Unidirectional laminate axial shear properties can be determined 
approximately from tensile tests of a [±45] unnotched laminate. 

In a [0 ./±45 ./90, ] laminate, on the other hand, the presence 
of many +45® fibers will significantly increase the shear stiff- 
ness of the laminate. Considerable tensile stresses will now be 
felt in the fiber direction of the ±45° layers. Also, due to the 
constraining effect of the ±45® layers, shear stress levels in the 
0 and 90® layers will be small and transverse normal stresses will 
dominate (see Appendix D) . Although layer axial tensile properties 
are less affected by fatigue loading than layer transverse tensile 
properties, the axial stress levels in the ±45® layers are also 
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considerably higher than the shear stresses in the 0 and 90® layers. 
Therefore, consideration must be given to the axial tensile and 
transverse tensile fatigue characteristics of layers as well as 
the axial shear fatigue characteristics. 

Example .- To illustrate the application of the foregoing 
fatigue model analysis to notched composite laminates, a boron/ 
epoxy l0^/±45] laminate will be analyzed. The laminate is assumed 
to contain a 2.54cm slit notch. Since at this time incomplete ex- 
perimental data are available, hypothetical variations of tensile 
and shear properties will be assumed. It will also be assumed in 
the present example that the axial Young's modulus of a layer, 
does not change significantly with cyclic loading. It will, how- 
ever, be assumed that the axial tensile strength of a unidirectional 
layer does exhibit some degradation as shown in Figure G-7. The 
curve shown for fatigue lifetime has been taken from experimental 
data presented in Reference 2, but the residual strength curves 
associated with the S-N curve are hypothetical. In axial shear, 

0 ^ 2 > behavior shown for residual strength and fatigue life- 

time has been approxmiated from tests presented in Reference 2 and 
from an S-N curve for epoxy material alone in tension shown in 
Reference 10. While an attempt was made to fabricate as realistic 
a curve as possible, these numbers must be regarded as purely 
hypothetical. The curve shown for variation of axial shear modulus 
G ^2 N has been estimated from scanty data presented in Refer- 

ence 2. Again, this curve must be regarded as hypothetical. 

The results of the failure analysis are shown in Figure G-8 
for the axial and transverse crack propagation modes. Increments 
for fatigue property variations were chosen as shown in Table G-1. 

It was assumed for simplicity’s sake that does not contribute to 
fatigue degradation. The results of the fatigue analysis for each 
layer are also shown in Table G-1. Predictions of laminate pro- 
perties from layer properties are presented along with resulting 
failure stresses in the axial and transverse crack propagation 
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modes. The variations of laminate axial shear properties with 
number of cycles are shown plotted in Figure G-8. 

For off-axis cracking, static failure analysis results for 
are shown in Figure G-9 for the laminate. Note that for off- 
axis analysis purposes, the laminate is now a [0/±452/90] laminate 
with a 1.9cm slit notch. To determine the degradation of x*y* 
shear properties during cyclic loading, the method of effective 
shear stress at the notch root was selected, equation (El). The 
value of the average effective notch shear stress as computed from 
equation (El) is shown in Table G-2 along with the results of the 
laminate and fatigue property degradation analyses for each incre- 
ment in cyclic loading considered. Also shown in this table are 
the resulting laminate properties at each incremental cycle, and 
the resulting residual strength of the notched composite in the 
off-axis propagation mode. Laminate properties as a function of 
number of cycles are also plotted in Figure G-9 to illustrate their 
effect on off-axis strength. It is noted that t®' , not t®* f is 
plotted in this figure, since it is the effective reduced shear 
strength which controls notched laminate behavior. 

As seen from the numbers tabulated in Tables G-1 and G-2, and 

the graphs in Figures G-8 and G-9, the off-axis propagation mode and 

the transverse crack propagation mode never occur for this laminate. 

The failure mode for initial static failure through the fatigue 

lifetime of the notched laminate is axial crack propagation. The 

resulting residual strength versus N curve is shown in Figure G-10. 

Note that, initially, the strength of the notched laminate increases 

with the number of cycles, due to a beneficial stress redistribution 

inside the laminate caused by decreasing shear modulus in 0® and 90® 

layers. Static residual failure in this region is due to a shear 

failure in the laminate which has been initiated by transverse 

2 

tensile failure in the -45® layers. After approximately 10 cycles, 
the residual strength of the laminate begins to decrease until it 
reaches the cyclic maximum stress at approximately 5 X 10^ cycles, 
which is the notched laminate fatigue lifetime, The decreasing 
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portion of the residual strength curve is due to laminate axial 
shear failure precipitated by axial tensile failure of the +45"^ 
layers . 

Figure G-11 shows how key elastic and strength properties de- 
grade for this particular laminate as a function of number of 
cycles. In this case, it is seen that the laminate axial modulus, 
and the axial shear modulus, do not vary significantly 

with fatigue loading. The maximum change is within a few percent. 

On the other hand, the axial shear strength has degraded consider- 
ably and is the main cause of laminate fatigue failure in this 
situation. Although this example is based upon hypothetical be- 
havior for layer fatigue, it serves to illustrate that a compre- 
hensive analysis can be performed provided that the proper input 
information is available. Also, the example shows that the 
assumption of a constant laminate axial Young's modulus remains 
valid as long as the axial Young's modulus of a layer does not 
significantly change with fatigue loading. 

As a means of assessing the effect of notch size on the 
static and fatigue failure behavior of a [ 02 /-^^] boron/epoxy lam- 
inate, the results of the static failure analysis for a 0.64cm 

notch with varying x® and G are presented in Figures G-12 and G-13. 

xy 

Although the fatigue example is not redone here, it is seen by 
comparing strength properties on a point by point laminate input 
property basis that the laminate with the smaller notch size is 
generally the stronger. A complete comparison for fatigue behavior 
must await valid input data and resulting re-analysis of both the 
2.54 and 0.64cm notch laminates. 
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TABLE G-1 


LAMINATE AND FATIGUE ANALYSIS RESULTS FOR 
t02/±45] B/Ep, 2.54 cm. NOTCH 

AXIAL/TRANSVERSE CRACK MODES 

S = 0. 85, = 243 MN/m2 



No. of 
Cycles 
N 

1 

102 

10 3 

10 

10® 

5 X 10® 

Laminate 

Elastic 

Properties 

GN/m^ 


112 

111.5 

110.4 

109.9 

109.3 

109.2 

°xy 

28.6 

28.4 


27.80 

27.60 

27.50 

Unnotched 

Laminate 

Failure 

Stresses 

MN/m2 

t 

a 

X 

710 

705 

697 

694 

689 

688 


285 

359 

351 

339 

284 

246 

Notched 

Laminate 

Failure 

Stresses 

MN/m2 

'^A 

250 

324 

317 

310 

255 

217 

'’SCFM 

545 

690 

676 

648 

538 

469 


Since <a^, trans crack never occurs 

SCFM X 
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TABLE G-2 


LAMINATE AND FATIGUE ANALYSIS RESULTS FOR 
[02/t45] B/Ep, 2.54 cm. NOTCH OFF-AXIS CRACK MODE 

S = 0.85, Tg' =161 MN/m2 



No. of Cycles 

1— I 

102 

10 

10^ 

10 8 

Laminate 

Elastic 

Properties 


65.55 

65.30 

64.43 

64.09 

63.45 

— 
Eyi (GN/m2) 

65.55 

65.30 


64.09 

63.45 

VI 1 1 

xy 

0.146 

0.144 

0.139 

0.137 

0.133 


18.88 

18.79 

18.44 

18.31 

18.05 

Unnotched 

Laminate 

Failure 

Stresses 

(MN/m2) 

266 

280 

315 

312 

306 

* T (MN/m2 ) 

416 

415 

411 

409 

397 

x' Axial 

Failure 

Stress 

a^, (MN/m2) 

175.8 

179.3 

184.8 

183.4 

179.3 

Laminate 

Tensile 

Failure 

Stress 

Off-Axis 

Mode 

0^5 (MN/m2) 

1 

351 .6 

358.5 

369.6 

366.8 

358.5 


*T® ' (not presented in this table) is used in failure 
computations. Figure G-9) . 
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AXIAL SHEAR STRENGTH 
r°j 10 MN/m^' 

AXIAL SHEAR MODULUS 
gn/m2 


ULTIMATE AXIAL SHEAR 
FAILURE STRAIN^ 

yjjl^,io-2 m/m 



NUMBER OF CYCLES, N 


FIGURE G-1. SHEAR PROPERTY DEGRADATION FOR 

HYPOTHETICAL MATERIALS A AND B. CYCLIC LOADING AT 
= 76 MN/M2 . 
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ULTIMATE AXIAL STRAIN^ ^ m/M 

(a) NOTCHLD LAMINATE TENSILE STRENGTH, (r . 
AXIAL CRACK FAILURE MODE. 



ULTIMATE AXIAL STRAIN. Vylt. M/M 

(b)COMPOSITE OVERSTRESS IN MATERIAL ADJACENT TO 
NOTCH. and corresponding FAR FIELD APPLIED 

GROSS STRESS, 

FIGURE G-2. STATIC AND FATIGUE BEHAVIOR OF BORON/EPOXY [01 
NOTCHED LAMINATE. NOTCH SIZE = 0.6ACM . = 400 MN/M^' 
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FIGURE G-5 RESIDUAL STRENGTH AND LIFETIME FOR 0.6A CM 

NOTCHED HYPOTHETICAL [0] B/Ep LAMINATE CYCLED AT S = 0.735. 
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AXIAL 

"crack" 

LENGTH 

a 


N = 103 102 1 



AXIAL 

"crack" 

LENGTH 

o 



APPLIED STRESS LEVEL <T , mn/m^ 


FIGURE 6-4. STATIC AXIAL "CRACK" GROWTH FOR 0.64 CM NOTCHED 
[01 B/EP LAMINATES COMPOSED OF HYPOTHETICAL MATERIALS A AND 
B HAVING PROPERTIES WHICH EXIST AFTER CYCLING hi o = 400 
MN/M^ FOR h 10^ AND 10^ CYCLES. 
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AXIAL 

"crack" 

LENGTH, 

ot 



FIGURE G>5. PREDICTED FATIGUE AXIAL "CRACK" GROWTH FOR 
0.6^1 CM NOTCHED [01 B/Ep LAMINATE. COMPOSED OF HYPOTHETICAL 
MATERIALS A AND B HAVING PROPERTIES WHICH EXIST AFTER 
CYCLING AT = AOO MN/M2. 


101 . 



LAMINATE 

AXIAL 

SHEAR 

STRENGTH 

mn/m^ 



(a) NOTCHED LAMINATE TENSILE STRENGTH, a 
AXIAL CRACK FAILURE MODE. ^ 


LAMINATE 

AXIAL 

SHEAR 

STRENGTH 

r° 

mn/m^ 



ULTIMATE AXIAL SHEAR STRAIN. Xjlt. M/M 


(b) COMPOSITE OVERSTRESS IN MATERIAL ADJ^^ 

TO NOTCH, and CORRESPONDING FAR FIELD 


APPLIED GROSS STRESS, 


FIGURE G-6. STATIC AND FATIGUE BEHAVIOR OF BORON/EPOXY 
[01 LAMINATE. NOTCH SIZE = 2.5A CM. 
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s 

0 r- 
9 
8 
7 

oL 



= 0 


NUMBER OF CYCLES^ N 

(a) AXIAL TENSILE STRENGTH, 




RGURE G-7. HYPOTHETICAL VARIATIONS OF FATIGUE BEHAVIOR 
AND LIFETIME OF (01 BORON/EPOXY 
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note: = 710 mn/m 
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FIfiURE 6-8. STATIC AND FATIGUE BEHAVIOR OF BORON/EPOXY IO 2 /+ 45 I NOTCHED LAMINATE. 

NOTCH SIZE = 2.5A CM 
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FIGURE G-9. OFF-AXIS FATIGUE BEHAVIOR OF |02/±45I SLIT NOTCHED LAMINATE. 

NOTCH SIZE = 2.54 cm 


OFF-AXIS 

RESIDUAL 

STRENGTH 

'^45 

mn/m^ 

AXIAL 

RESIDUAL 

STRENGTH 

mn/m^ 


FIGURE 6-10. 



NUMBER OF CYCLES^ N 

residual strength of notched [02/±451 BORON/ 

EPOXY U\F1INATE CYCLED AT S = 0.85 . 

NOTCH SIZE = 2.54 cm 
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400 _ 
300 _ 
200 _ 
100 _ 



PREDICTED LAMINATE AXIAL MODULI AND SHEAR 
STRENGTH OF NOTCHED I02/±A51 BORON/EPOXY 

LAMINATE CYCLED AT S = 0.85. 

NOTCH SIZE = 2.5A CM 


FIGURE G-11. 
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XV LAMINATE SHEAR MODULUS 

6^1 yi (gn/m^) 

FIGURE 6-13. OFF-AXIS FATIGUE BEHAVIOR OF BORON/EPOXY 
102/±A51 LAMINATE with slit notch. NOTCH SIZE = 0.6A 

(EQUIVALENT NOTCH IN X'Y‘ COORDINATE SYSTEM = 0.A5 CM) 
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LIBRARY-CARD ABSTRACT 


This report describes the development of 
a semi-empirical analysis for prediction and 
correlation of fatigue crack growth, residual 
strength, and fatigue lifetime for fiber com- 
posite laminates containing notches. The analy- 
sis is based upon experimentally observed failure 
modes and predicts notched laminate failure from 
unidirectional layer tensile and shear fatigue 
properties. Examples are given for [0] and multi- 
directional laminates which illustrate applica- 
tion of the method. 



